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PREFACE 


The  work  reported  herein  was  conducted  as  part  of  the  Water  Quality 
Research  Program  (WQRP) .  The  WQRP  Is  sponsored  by  the  Headquarters ,  US  Army 
Corps  of  Engineers  (HQUSACE) ,  and  Is  assigned  to  the  US  Army  Engineer 
Waterways  Experiment  Station  (WES)  under  the  purview  of  the  Environmental 
Laboratory  (EL) .  Funding  was  provided  under  Department  of  the  Army  Appropria¬ 
tion  96X3121,  General  Investigation.  The  WQRP  Is  managed  under  the  Environ¬ 
mental  Resources  Research  and  Assistance  Programs  (ERRAP) ,  Mr.  J.  L.  Decell, 
Manager.  Mr.  Robert  C.  Gunkel  was  Assistant  Manager,  ERRAP,  for  the  WQRP. 
Technical  Monitors  during  this  study  were  Mr.  David  Buelow,  Mr.  Jim  Gottesman, 
and  Dr.  John  Bushman,  HQUSACE. 

This  report  was  prepared  by  the  Hydraulics  Laboratory  (HL) ,  WES,  and  the 
St.  Anthony  Falls  Hydraulic  Laboratory  (SAFHL)  of  the  Department  of  Civil  and 
Mineral  Engineering,  University  of  Minnesota,  Minneapolis,  as  a  summary  of 
results  from  research  Into  the  area  of  gas  transfer  at  low-head  hydraulic 
structures.  This  effort  was  funded  under  WQRP  Work  Unit  32369  entitled 
"Reaeratlon  at  Low-Head  Structures.”  The  study  was  conducted  under  the 
general  direction  of  Messrs.  Frank  A.  Herrmann,  Jr.,  Director,  HL;  Richard  A. 
Sager,  Assistant  Director,  HL;  and  Glenn  A.  Pickering,  Chief,  Hydraulic 
Structures  Division  (HSD) ,  HL. 

This  report  was  prepared  by  Mr.  Steven  C.  Wilhelms  of  the  Reservoir 
Water  Quality  Branch  (RWQB) ,  HSD;  Dr.  John  S.  Gulliver,  SAFHL,  and  Mr.  Kenneth 
Parkhlll,  SAFHL,  under  the  direct  supervision  of  Dr.  Jeffery  P.  Holland,  for¬ 
mer  Chief,  RWQB.  Mr.  Perry  Johnson,  US  Bureau  of  Reclamation,  and  Dr.  Charles 
Bohac,  Tennessee  Valley  Authority,  provided  review  comments.  Field  and  liter¬ 
ature  studies  were  conducted  by  personnel  from  HL  and  SAFHL.  Assistance  In 
data  analysis  and  report  preparation  was  provided  by  Ms.  Laurln  I.  Yates, 

RWQB,  and  Mrs.  Barbara  A.  Parsons,  HSD,  and  Mr.  Benjamin  Erickson,  SAFHL.  The 
report  was  edited  by  Mrs.  Marsha  C.  Gay,  Information  Technology  Laboratory, 
WES. 

At  the  time  of  publication  of  this  report.  Director  of  WES  was 
Dr.  Robert  W.  Whalln.  Commander  and  Deputy  Director  was  COL  Leonard  G. 
Hassell,  EN.  Dr.  Roger  L.  A.  Arndt  was  Director,  SAFHL. 

This  report  should  be  cited  as  follows: 
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CONVERSION  FACTORS,  NON-SI  TO  SI  (METRIC) 
UNITS  OF  MEASUREMENT 


Non-SI  units  of  measurement  used  in  this  report  can  be  converted  to 
(metric)  unit  as  follows; 

Multiply  Bv 

degrees  (angle)  0.017AS329 

feet  0.3048 


SI 

To  Obtain 

radians 

metres 
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REAERATION  AT  LOW-HEAD  HYDRAULIC  STRUCTURES 


PART  I :  INTRODUCTION 

Background 

1.  Presently,  one  of  the  most^  cited  water  quality  parameters  in  the 
freshwater  hydrosphere  (rivers,  lakes,  and  reservoirs)  is  dissolved  oxygen 
(DO) .  The  oxygen  concentration  in  surface  waters  is  a  prime  indicator  of  the 
quality  of  that  water  for  human  use  as  well  as  use  by  the  aquatic  biota.  Many 
naturally  occurring  biological  and  chemical  processes  use  oxygen,  thereby 
diminishing  the  DO  concentration  in  the  water.  The  physical  process  of  oxygen 
transfer  or  oxygen  absorption  from  the  atmosphere  or  air  bubbles  acts  to  re¬ 
plenish  the  used  oxygen.  This  process  is  termed  reaeration. 

2.  Low-head  hydraulic  structures  within  the  US  Army  Corps  of  Engineers 
are  generally  associated  with  navigation  projects.  These  structures  are  usu¬ 
ally  "run-of-the-river"  and  are  used  to  maintain  a  constant  upstream  pool 
elevation.  The  oxygen  transfer  in  these  deeper,  slower  pools  is  lower  than 
that  of  the  open  river.  Biological  and  chemical  oxygen  demands  may  accumulate 
and  concentrate  in  the  impoundment  and  thereby  degrade  the  DO  concentration  in 
the  stored  water  because  of  the  excess  demand  compared  to  reaeration  poten¬ 
tial.  Without  sufficient  reaeration,  release  of  this  water  may  pose  an 
environmental  and  water  quality  concern. 

3.  Some  hydraulic  structures  exhibit  remarkable  reaeration,  while 
others  do  very  little  to  increase  DO.  If  the  DO  in  releases  is  lower  than 
desired  or  required,  operational  or  structural  modifications  or  artificial 
aeration  can  be  employed  to  improve  the  release  DO  concentration.  The  design 
engineer  will  be  faced  with  the  need  to  evaluate  the  oxygen  transfer  charac¬ 
teristics  of  existing  conditions  at  a  hydraulic  structure  for  comparison  with 
the  characteristics  of  a  proposed  modification.  As  in  any  engineering  appli¬ 
cation,  it  is  imperative  to  possess  a  clear  definition  of  the  processes  that 
affect  water  quality  and  the  expected  change  as  a  result  of  implementing  the 
particular  improvement  technique.  Hence,  the  potential  impacts  of  alternative 
techniques  on  the  reaeration  processes  must  be  thoroughly  understood  and 
quantified. 

4.  In  the  past,  the  focus  of  interest  in  gas  exchange  at  hydraulic 
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structures  has  been  the  transfer  of  atmospheric  gases:  absorption  of  oxygen 
to  replace  DO  used  by  aquatic  processes  and  the  absorption  of  nitrogen, 
potentially  resulting  In  nitrogen  supersaturation.  More  recently,  however, 
the  desorption  of  volatile  organics  or  toxics  that  may  be  dissolved  In  the 
water  has  become  Important.  The  air-water  transfer  of  any  chemical  Is  called 
gas  transfer  because  the  chemical  Is  a  gas  In  one  of  the  two  phases.  In  gen¬ 
eral,  the  physical  processes  that  Influence  oxygen  absorption  also  affect  the 
transfer  of  any  dissolved  volatile  compound.  Thus,  It  becomes  even  more 
Important  for  the  design  engineer  to  possess  a  thorough  understanding  of 
(a)  the  physics  of  gas  transfer,  (b)  the  quantification  of  gas  transfer, 

(c)  the  Important  physical  processes,  and  (d)  the  hydraulic  conditions  that 
can  enhance  or  degrade  gas  transfer. 

Objective  and  Scope 

5.  The  objective  of  this  report  Is  to  familiarize  field  engineers  with 
the  oxygen  transfer  process  and  the  oxygen  transfer  characteristics  of  various 
low-head  hydraulic  structures.  The  physics  of  gas  transfer  are  conceptually 
explained  and  the  mathematical  description  of  the  gas  transfer  process  Is 
developed.  The  Important  physical  processes  and  their  Impact  on  the  variables 
In  the  gas  transfer  equation  are  Identified.  The  hydraulic  conditions  that 
contribute  to  these  physical  processes  are  described  and  applied  to  oxygen 
transfer,  or  reaeratlon. 

6.  It  Is  hoped  that  field  engineers,  when  familiar  with  the  conceptual 
descriptions  provided  in  this  report,  can  qualitatively  evaluate  the  oxygen 
transfer  characteristics  at  a  structure  based  solely  on  observed  hydraulic 
conditions,  e.g.,  they  will  be  able  to  estimate  whether  the  structure,  upon 
testing,  would  exhibit  a  large  or  small  degree  of  gas  transfer.  With  the 
mathematical  description  of  oxygen  transfer  at  a  "generic"  type  of  structure 
and  an  understanding  of  how  hydraulic  conditions  contribute  to  oxygen  trans¬ 
fer,  field  engineers  should  be  able  to  "bracket"  or  roughly  estimate  the 
transfer  that  is  occurring  at  a  specific  structure. 

7.  The  hydraulic  structures  at  most  low-head  projects  usually  consist 
of  a  gated  sill,  gated  low-head  spillway,  and  a  fixed-  or  adjustable-crest 
weir.  The  gas  transfer  analyses  reported  herein  emphasize  these  "generic" 
types  of  structures.  More  complete  descriptions  of  the  geometry  and  flow 
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conditions  at  these  structures  are  given  in  Part  IV,  "Hydraulics  at  Various 
Structv^res . " 
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PART  II:  REAERATION  THEORY 


Physics  of  Gas  Transfer 


8.  The  transfer  of  any  chemical  across  an  air-water  interface  can  be 
considered  two  processes  acting  together:  (a)  transport  due  to  molecular 
diffusion  and  (b)  transport  due  to  turbulent  mixing.  These  are  characterized 
by  the  capture  (or  release)  of  gas  molecules  at  the  air-water  Interface  and 
subsequent  distribution  of  the  gaa  throughout  the  water  body.  Molecular  dif¬ 
fusion  is  caused  by  the  inherent  kinetic  energy  possessed  by  the  gas  mole¬ 
cules,  whereas  transport  due  to  turbulent  mixing  results  from  the  application 
of  external  forces  on  the  water  body. 


Figure  1.  Transport  due  to 
diffusion  in  quiescent 
water  (after  Tsivoglou  and 
Wallace  1972) 


9.  Consider  the  quiescent  body  of 
water  shown  in  Figure  1.  If  a  group  of 
dissolved  molecules  could  be  placed  in  the 
water  body  without  disturbing  the  water, 
the  dissolved  molecules  would  gradually 
spread  throughout  the  water  body  and  even¬ 
tually  achieve  a  uniform  distribution.  The 
molecules  would  accomplish  this  in  totally 
still  water  because  of  their  inherent 
kinetic  energy  associated  with  their  sur¬ 
rounding  temperature.  The  movement  of  the 
molecules,  which  is  entirely  random,  causes 
the  molecules  to  uniformly  distribute 
throughout  the  water  body.  This  process  is 


termed  molecular  diffusion. 

10.  Molecular  diffusion  is  usually  described  by  Fick's  first  law* 


d  =  -D„ 

ID 


dC 

Hi 


(1) 


which  states  that  J  ,  the  mass  flux  per  unit  area  across  an  interface  or  a 
plane  (in  a  direction  away  from  the  higher  concentration,  hence,  the  negative 


*  For  convenience,  symbols  and  unusual  abbreviations  are  listed  and  defined 
in  the  Notation,  Appendix  E. 
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sign) ,  is  proportional  to  the  concentration  gradient  dC/ds  across  the 
Interface,  where  dC  Is  the  change  In  concentration  across  the  Interface 
thickness  ds  .  D„  Is  the  molecular  diffusion  coefficient  and  Is  a  measure 
of  the  ability  of  a  medium  to  diffuse  dissolved  material.  This  relationship 
Implies  that  even  though  the  molecular  diffusion  coefficient  may  be  small,  the 
flux  from  one  volume  with  a  high  dissolved  concentration  into  a  volume  with 
zero  dissolved  concentration  will  initially  be  very  high  because  of  the  large 
gradient  across  the  Interface  between  the  volumes. 

11.  For  an  example  application  of  these  concepts  to  gas  transfer,  again 
consider  a  completely  quiescent  homogeneous  body  of  water.  Assume  that  the 
body  of  water  Is  devoid  of  oxygen  and  that  the  surface  Is  Instantaneously 
exposed  to  the  atmosphere.  Initially,  there  would  be  a  very  high  rate  of 
oxygen  absorption  into  a  very  thin  layer  at  the  water  surface  because  of  the 
large  concentration  gradient  across  the  Interface  (Pick's  law).  However, 
since  this  Is  a  quiescent  water  body,  only  molecular  diffusion  will  cause 
oxygen  molecules  to  move  down  into  the  water  body.  Since  the  diffusion  down¬ 
ward  is  impeded  by  the  molecular  structure  and  viscosity  of  the  water,  oxygen 
molecules  will  tend  to  collect  rapidly  in  the  surface  layers.  Consequently, 
the  rate  of  absorption  would  fall  quickly  because  of  the  rapid  decrease  in  the 
concentration  gradient  across  the  air-water  Interface .  Thus ,  the  rate  of  gas 


or  mlxlnt 


12.  Consider  that  same  body  of 


water,  devoid  of  oxygen  and  with  its 
surface  instantaneously  exposed  to  the 
atmosphere.  However,  this  time,  the  water 
is  being  agitated  causing  turbulent  mixing 
(Figure  2) .  Turbulent  mixing  transports  a 
water  element  to  the  water  surface,  which 
degrades  the  thin  film  at  the  surface  for  a 
short  period.  Although  the  period  of  expo¬ 
sure  to  the  atmosphere  may  be  very  short 
and  the  molecular  diffusion  coefficient  may 
be  small,  the  element  will  absorb  a  signif¬ 
icant  amount  of  oxygen  because  of  the  large 
gradient  across  the  interface  (according  to 


Figure  2.  Transport  due  to 
diffusion  and  turbulent 
mixing 
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Pick's  law,  the  rate  of  transfer  would  be  large).  Similarly,  when  that 
element  is  transported  by  turbulent  mixing  into  the  body  of  fluid,  the  rate  of 
oxygen  transfer  to  surrounding  elements  will  be  very  high  because  of  the 
relatively  large  gradient.  It  is  easily  deduced  from  this  description  that 
turbulent  mixing  greatly  speeds  oxygen  transfer. 

13.  It  is  Important  to  realize  that  gas  transfer  is  governed  by  the 
same  physical  processes,  turbulent  mixing  and  molecular  diffusion,  regardless 
of  the  direction  of  transfer,  i.e.,  gas  absorption  or  desorption.  This  fact 
provided  the  basis  for  developing  all  tracer  gas  technologies.  Tsivoglou 
et  al.  (1965)  and  Tsivoglou  et  al.  (1968)  reported  this  basis  and  the  develop¬ 
ment  of  a  technique  that  used  radioactive  krypton-85,  injected  into  the  water, 
as  a  tracer  for  oxygen  absorption.  The  krypton-85  was  being  transferred 
(desorbed)  to  the  atmosphere,  while  oxygen  was  being  absorbed.  Rathbun  et  al. 
(1978)  modified  the  radioactive  tracer  technique  with  the  injection  of  hydro¬ 
carbon  gases,  such  as  propane  and  ethylene,  as  tracers.  Uilcock  (1983)  used 
methyl  chloride  as  a  tracer  gas  in  this  same  fashion.  McDonald,  Gulliver,  and 
Wilhelms  (1990)  reported  preliminary  measurements  of  gas  transfer  with  in  situ 
methane  and  sulfur  hexafluoride.  McDonald  and  Gulliver  (1991)  presented  the 
development  of  methane  as  an  alternative  gas  for  measuring  the  gas  transfer  in 
a  particular  flow  situation  at  hydraulic  structures.  The  critical  point  upon 
which  all  this  work  was  based  is  that  measurements  made  with  one  gas  can  be 
used  to  calculate  the  transfer  of  another  gas  (Gulliver,  Thene,  and  Rindels 
1990). 


Mathematical  Description 

14.  Gas  transfer  is  usually  considered  to  be  a  first-order  process  in 
which  the  rate  of  change  of  the  chemical's  concentration  in  the  water  is 
linearly  dependent  on  the  ambient  concentration.  The  driving  force  in  the 
transfer  process  is  the  difference  between  the  actual  concentration  of  the 
dissolved  gas  in  the  water  and  the  concentration  that  corresponds  to  equilib¬ 
rium  with  the  air.  This  equilibrium  concentration  or  "saturation  concentra¬ 
tion"  is  defined  using  a  Henry's  law  constant  as  follows: 

=  HiPi  (2) 
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where 


Cal  **  saturation  concentration  of  gas  1  ,  g/m’ 

Hi  •  Henry's  law  constant  for  gas  1  ,  g/a^  -  atn 

?!  -  partial  pressure  of  gas  1  In  the  atmosphere,  atm 

Henry's  law  states  that  at  a  given  temperature,  a  liquid  can  absorb  an  amount 
of  gas  that  Is  proportional  to  the  partial  pressure  of  that  gas  In  the  overly¬ 
ing  atmosphere.  Thus,  an  "equilibrated"  state  exists  at  the  saturation  con¬ 
centration.  There  Is  a  rate  coefficient  on  the  air  side  (gas  film  coeffi¬ 
cient)  and  on  the  water  side  (liquid  film  coefficient)  of  the  Interface.  The 
Inverse  of  these  coefficients  can  be  seen  as  "resistance"  to  gas  transfer 
(Llss  and  Slater  1974) ,  where  the  air-  and  water-side  resistance  act  In 
series.  For  most  compounds,  such  as  O2.  CO2,  CH4,  N2,  etc.,  the  water-side 
resistance  Is  much  greater  than  the  alr-slde  resistance,  and  only  the  liquid 
film  coefficient  Is  considered.  The  gas  film  coefficient  can  also  be  Impor¬ 
tant  in  the  transfer  of  H2O  (evaporation),  NH3,  SO2,  and  some  herbicides  and 
pesticides  designed  for  this  purpose.  Because  this  report  on  reaeratlon  con¬ 
cerns  oxygen  transfer,  only  the  water-side  resistance  need  be  considered. 

15.  A  measure  of  the  driving  force  causing  oxygen  transfer  can  there¬ 
fore  be  defined  as  the  difference  between  the  concentration  in  the  water  and 
the  saturation  concentration.  This  quantity  Is  called  the  "saturation 
deficit"  and  is  mathematically  defined  by 


D  -  C,  -  C 


(3) 


where 

D  -  saturation  deficit 
C,  —  saturation  concentration 
C  -  ambient  concentration  In  water 

At  saturated  conditions,  the  deficit  is  zero.  If  the  deficit  Is  positive,  the 
water  Is  "undersaturated."  If  the  deficit  Is  negative,  then  the  water  Is 
"supersaturated."  Oxygen  would  be  absorbed  for  the  former  and  desorbed  for 
the  latter  to  achieve  the  equilibrated  saturated  state. 

16.  To  develop  an  expression  that  describes  oxygen  transfer  as  a 
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first-order  process,  the  flux  across  the  air-water  interface  is  typically 
written  as 


J  -  kL(C,  -  C)  =  ki.  D 


(4) 


where  kj,  is  the  liquid  film  coefficient  for  oxygen  transfer.  Examination  of 
this  equation  and  Equation  1  indicates  that  the  liquid  film  coefficient  is 
related  to  the  concentration  gradient  near  the  air-water  interface,  and  may  be 
Interpreted  as  the  ease  with  which  dissolved  gases  move  across  the  air-water 
interface.  The  liquid  film  coefficient  is  dependent  on  the  internal  structur 
or  ordering  of  water  molecules  and  the  breakdown  of  that  structure  by  mixing. 

17 .  To  determine  the  rate  of  change  in  concentration  for  a  well-mixed 
volume  of  water,  the  total  rate  of  mass  flux  into  a  water  body  must  be  calcu¬ 
lated  by  multiplying  Equation  4  by  the  surface  area  and  then  dividing  this 
quantity  by  the  volvune  of  water.  Dividing  the  total  mass  flux  by  the  volume 
results  in  the  rate  of  change  in  concentration: 


JA  .  dC 

T  Bt 


dD 


=  kL^(C.  -  C) 


(5) 


when  C,  is  constant  and  where 

A  -  surface  area  associated  with  the  volume  V  ,  over  which  transfer 
occurs 

V  -  volxime  of  the  water  body  over  which  A  is  measured 
dC/dt  -  rate  of  change  of  concentration 
dD/dt  -  rate  of  change  of  the  saturation  deficit 

18.  Assuming  that  kj,  ,  A  ,  and  V  are  constant  over  the  time  of 
flow.  Equation  5  can  be  integrated  to  the  following  mathematical  model  of 
water-side  controlled  gas  transfer: 


Df 

T5I 


exp  (-k^  ^  t) 


=  exp  (-kLat)  =  exp  (-Kjt) 


(6) 
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where 


D;,D^  -  final  and  initial  oxygen  deficits,  respectively 
t  -  elapsed  time  from  initial  to  final  deficits 
a  -  specific  surface  area,  A/V 
K2  -  reaeration  or  oxygen  transfer  coefficient 

From  the  perspective  of  flow  in  a  stream  reach  or  through  a  hydraulic  struc¬ 
ture,  the  initial  and  final  deficits  would  be  the  upstream  and  downstream 
deficits,  respectively,  and  the  elapsed  time  would  be  the  time  of  flow  from 
the  upstream  to  downstream  locations.  As  presented  later,  some  researchers 
have  defined  a  quantity  r  as  the  deficit  ratio,  which  is  the  Inverse  of 
Equation  6  such  that 


r 


exp  (  kj^at  ) 


(7) 


However,  the  mathematical  qualities  of  this  formulation  are  less  desirable 
than  those  of  Equation  6  (Wilhelms  and  Smith  1981;  Rlndels  1990). 

19.  A  convenient  parameter  called  transfer  efficiency  can  be  defined  as 
the  fractional  part  of  the  incoming  deficit  that  is  satisfied  as  the  water 
flows  through  the  structure  or  stream  reach.  This  efficiency  E  can  be 
expressed  mathematically  with  Equation  6  as 


E 


(G^  -  Ci) 

“CO 


1  -  exp  (-kj,  ^  t) 


(8) 


E  *  1  -  exp  (-Kjt) 


(9) 


where  Cf  and  Cj  are  the  final  and  initial  oxygen  concentrations,  respec¬ 
tively.  If  the  transfer  efficiency  is  zero,  then  there  was  no  oxygen  transfer 
and  the  downstream  concentration  equals  the  upstream  concentration.  If  the 
transfer  efficiency  equals  1.0,  then  all  of  the  upstream  deficit  was  met  by 
the  oxygen  transfer  and  the  downstream  concentration  is  at  saturation.  Fur¬ 
thermore,  the  transfer  efficiency  of  oxygen  can  be  transferred  to  other  gases 
and  vice  versa  through  the  indexing  method  developed  by  Gulliver,  Thene,  and 
Rindels  (1990). 
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Important  Physical  Processes 


20.  Three  processes  occur  at  hydraulic  structures  that  can  signifi¬ 
cantly  increase  oxygen  transfer.  They  can  be  related  to  the  parameters  given 
in  Equations  8  and  9.  The  impacts  of  these  processes  are  governed  mainly  by 
the  fluid  mechanics  and  flow  conditions  at  the  structure.  The  following  pro¬ 
vides  a  general  description  of  these  processes  and  their  effect  on  gas 
transfer  at  hydraulic  structures: 

a.  Turbulent  mixing  at  the  water  surface  and  within  the  body  of  the 
flowing  water.  It  would  seem  logical  that  the  rate  of  turbulent 
mixing  would  significantly  affect  gas  transfer  because  of  the 
concept  of  water-surface  renewal  (water  surface  that  is  swept 
away  from  the  surface  and  "renewed"  with  water  from  below) 
(Danckwerts  1951)  causing  increased  gas  transfer.  A  high  degree 
of  turbulent  mixing,  such  as  occurs  on  the  face  of  a  spillway  or 
in  a  tailwater  plunge  pool,  would  increase  Ul  and  likewise 

E  . 

b .  Increased  interfacial  area  resulting  from  air  that  has  been 
entrained  into  the  flow.  When  air  is  entrained  into  the  flow 
either  from  the  surface  or  at  a  plunge  point,  the  surface  area 
available  for  gas  transfer  can  increase  dramatically.  Gulliver, 
Thene,  and  Rindels  (1990)  estimated  that  entrained  air  due  to 
free  surface  aeration  in  a  l-ft*-deep  flow  on  a  30-degree  slope 
can  increase  the  air-water  surface  area  by  a  factor  of  nearly 
500  compared  to  the  unit  area  of  surface  exposed  to  the  atmo¬ 
sphere.  Thus,  if  air  is  entrained,  gas  transfer  should  increase 
significantly  for  a  given  flow  condition.  Although  not  defined, 
the  amount  of  entrained  air  should  also  be  a  factor  in  estab¬ 
lishing  the  oxygen  transfer  character  of  a  hydraulic  structure. 
The  water-surface  area  A  is  greatly  increased  when  air  is 
entrained  into  the  flow.  This  would  also  act  to  increase  E  . 

c .  Increased  saturation  concentration  from  the  higher  pressure 
experienced  bv  bubbles  in  the  plunpe  pool.  In  addition  to  the 
contribution  that  air  bubbles  make  to  the  air-water  surface 
area,  absorption  of  atmospheric  gases  from  the  air  bubbles  can 
be  increased  because  of  the  increased  pressure  that  the  bubbles 
experience  as  they  are  transported  into  the  depth  of  the 
structure's  stilling  basin.  Increased  hydrostatic  pressure  on 
entrained  air  causes  an  increase  in  the  saturation  concentration 
(Equation  2)  and  thereby  increases  the  saturation  deficit 
(Buck,  Miller,  and  Sheppard  1980;  Wilhelms,  Schneider,  and 
Howington  1987;  Wilhelms  and  Gulliver  1990;  McDonald  1990). 

Thus,  although  the  transfer  efficiency  does  not  change  due  to 
this  effect,  the  downstream  concentration  would  increase  because 
of  an  increase  in  . 


*  A  table  of  factors  for  converting  non-SI  units  of  measurement  to  SI 
(metric)  units  is  found  on  page  5. 
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Each  of  these  processes  is  Included  either  directly  or  indirectly  in  Equa¬ 
tions  8  and  9.  Turbulent  mixing  is  characterized  by  the  liquid  film  coeffi¬ 
cient  ki,  .  The  increase  in  oxygen  transfer  due  to  greater  interfacial  area 
because  of  entrained  air  bubbles  is  Included  in  the  interfacial  area  term  A 
The  effects  of  pressure  on  oxygen  transfer  in  a  plunge  pool  result  in  an  in¬ 
creased  saturation  concentration  C,  .  This  Increase  in  C,  can  cause  super¬ 
saturated  oxygen  and  nitrogen  concentrations,  compared  to  the  surface  satura¬ 
tion  concentration.  By  recognizing  these  processes  and  their  impacts  on 
oxygen  transfer,  flow  conditions  observed  in  a  physical  model  or  full-scale 
project  can,  at  a  very  minimum,  be  qualitatively  evaluated.  If  one  objective 
in  the  operation  of  a  hydraulic  structure  is  to  provide  for  oxygen  transfer, 
this  evaluation  can  identify  means  of  meeting  that  objective. 
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PART  III:  MEASUREMENT  TECHNIQUES 


21.  The  efficiency,  along  with  the  variables  in  the  exponents  of  Equa¬ 
tions  8  and  9,  defines  the  gas  transfer  characteristics  of  a  structure.  To 
calculate  efficiency  or  estimate  those  variables,  upstream  and  downstrefim 
dissolved  gas  concentrations  are  required.  Essentially,  two  methods  are 
available  for  collection  of  these  data:  (a)  direct  in  situ  measurements  and 
(b)  discrete  sample  collection  with  a  subsequent  analysis  for  gas  concentra¬ 
tions.  These  two  methods  are  discussed  in  the  following  paragraphs  and  an 
uncertainty  analysis  is  presented  that  defines  the  minimum  values  required  for 
an  accurate  assessment.  Additionally,  conditions  for  testing  a  hydraulic 
structure  are  discussed  that  will  result  in  a  complete  picture  of  the  oxygen 
transfer  character  of  the  structure. 

In  Situ  Measurements 

22.  In  situ  measurements  are  usually  made  for  DO  with  a  portable 
electronic  meter  that  uses  a  polarographic  probe  (Yellow  Springs  Instrument 
Company  1975;  American  Public  Health  Association  1975).  The  probe  consists  of 
two  metal  electrodes  in  contact  with  an  electrolyte  separated  from  the  test 
water  by  a  membrane  that  is  permeable  to  oxygen.  A  polarizing  voltage  applied 
across  the  electrodes  causes  oxygen  that  has  permeated  the  membrane  to  react 
at  the  cathode,  resulting  in  an  electrical  current  flow.  Oxygen  passes 
through  the  membrane  at  a  rate  proportional  to  its  concentration  in  the  water. 
If  the  concentration  is  high,  more  oxygen  passes  through  the  membrane,  result¬ 
ing  in  an  increased  current  flow.  The  current  flow  is  calibrated  to  corre¬ 
spond  to  a  specific  DO  concentration. 

23.  These  probes  are  usually  very  reliable,  simple  to  operate  and  main¬ 
tain,  and  relatively  easy  to  calibrate.  Failure  to  provide  adequate  measure¬ 
ments  has  occurred  only  under  physical  or  biological  fouling  of  the  probe  or 
extremely  cold  weather  conditions.  They  are  accurate  to  1  percent  of  full 
scale  or  0.1  mg/i.  They  may  be  calibrated  with  an  air  contact  or  chemical 
titrametric  method.  Typically,  these  instruments  incorporate  a  thermistor  for 
temperature  measurements  and  a  stirrer  for  moving  water  past  the  permeable 
membrane . 
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Discrete  Sampling  and  Analysis  Techniques 


24.  Discrete  sample  collection  with  subsequent  analysis  Is  usually  loh- 
plemented  when  gases  such  as  H2S  or  CO2  are  of  Interest,  although  DO  can  also 
be  determined  with  a  discrete  sampling  method.  A  discrete  sampling  technique 
may  also  be  used  when  oxygen  concentrations  are  at  saturation  and  a  tracer  gas 
must  be  used  to  determine  the  transfer  characteristics.  There  are  several 
methods  for  collection  and  analysis  of  discrete  samples.  The  method  used 
depends  upon  the  dissolved  gas  of  Interest.  A  chemical  tltrametrlc  technique 
can  be  used  for  accurate  DO  measurements.  Gas  chromatography  can  be  used  for 
DO,  nitrogen,  hydrogen  sulfide,  carbon  dioxide,  methane,  and  any  other  vola¬ 
tile  gas  of  sufficient  concentration.  Chemical  tltrametrlc  and  colorimetric 
methods  can  also  be  used  for  some  of  these  gases,  with  variable  measurement 
accuracy . 

Sampling  techniques 

25.  The  sampling  instrument  usually  depends  upon  the  gas  of  Interest 
and  the  analysis  technique.  For  example,  samples  for  tltrametrlc  analysis  for 
DO  are  usually  collected  in  300-ml  biochemical  oxygen  demand  (BOD)  bottles  in 
a  sampler  specifically  designed  for  this  purpose.  Thene  and  Gulliver  (1989) 
designed  a  similar  sampler  for  use  with  AO-ml  bottles  In  propane  and  methane 
tracer  gas  analyses.  For  any  sampling  technique,  the  critical  Issue  Is  to 
collect  and  analyze  a  sample  that  Is  "undisturbed."  relative  to  in  situ  mea¬ 
surements.  Thus,  care  must  be  exercised  In  the  collection  and  handling  of 
samples  to  ensure  that  no  loss  or  gain  of  gas  occurs  as  a  result  of  the 
sampling  process.  This  may  dictate  that  a  special  sampler  be  designed  and 
constructed  to  meet  the  needs  of  the  particular  study. 

Sample  analysis  techniques 

26.  As  previously  mentioned,  a  chemical  tltrametrlc  technique  is  avail¬ 
able  for  determining  the  DO  concentration.  For  the  azide  modification  of  the 
Winkler  titration  technique,  divalent  manganese  sulfate  and  alkaline  lodide- 
azlde  are  added  to  the  300-ml  BOD  bottles.  They  chemically  react,  producing  a 
divalent  manganous  hydroxide  precipitate.  The  oxygen  In  the  water  sample 
rapidly  oxidizes  an  equivalent  amount  of  the  manganous  hydroxide  to  hydroxides 
of  higher  valency  states.  Adding  sulfuric  acid  to  this  solution  causes  the 
reversion  of  the  oxidized  manganese  to  the  divalent  state,  while  liberating 
iodine  equal  to  the  original  DO  content.  The  DO  is  determined  by  titrating 
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the  iodine  in  the  solution  to  Its  end  point  with  standard  solutions  of  sodium 
thiosulfate.  More  detail  is  given  in  the  handbook  Standard  Methods  for  the 
Examination  of  Water  and  Wastewater  (American  Public  Health  Association  1975) . 

27.  This  technique  is  relatively  simple,  but  is  more  time-consuming 
than  the  polarographic  probe  technique.  Care  must  be  exercised  in  sample 
handling  to  prevent  gas  transfer,  but  in  general,  this  technique  results  in 
measurements  that  are  reliable  and  accurate. 

28.  Gas  chromatography  is  an  analysis  method  useful  for  DO,  nitrogen, 
hydrogen  sulfide,  carbon  dioxide,  propane,  or  methane.  The  technique  is  based 
on  "sieving"  component  gases  in  a  sample  with  a  "molecular  sieve"  or  absorbent 
polymer.  The  gas  sample  passes  through  a  long  tube  (column)  packed  with  a 
porous  solid.  The  solid  acts  like  a  sieve,  slowing  the  movement  of  larger 
diameter  gas  molecules  or  absorbing  particular  gases.  This  results  in  separa¬ 
tion  of  the  sample  component  gases  of  interest  by  the  time  the  sample  exits 
the  column  at  the  detector.  Output  from  the  detector  is  calibrated  to  corre¬ 
spond  to  a  certain  amount  of  gas .  The  t3rpe  of  detector  used  in  the  gas 
chromatograph  (GC)  depends  upon  the  gas  of  interest.  For  DO  and  nitrogen,  a 
thermal  conductivity  detector  can  be  used  in  the  GC.  For  propane  and  methane, 
a  flame  ionization  detector  is  used  in  the  GC.  For  sulfur  hexafluoride,  an 
electron  capture  device  is  used  in  the  GC.  These  are  examples  of  the  sampled 
gases  and  tracer  gases  of  interest  for  hydraulic  structures. 

29.  In  using  the  GC  to  determine  the  mass  of  a  volatile  compound,  a 
headspace  sample  preparation  technique  has  proven  successful  (Thene  1988).  A 
headspace  of  inert  gas  is  created  in  the  sample  bottle.  By  agitating  the 
sample  bottle,  nearly  all  the  gas  in  the  water  is  stripped  into  the  headspace. 
Samples  from  this  headspace  volume  are  extracted  and  injected  into  the  GC  for 
determination  of  compound  mass.  The  actual  gas  concentration  in  the  water 
sample  is  calculated  based  on  the  mass  of  the  compound  in  the  headspace,  the 
relative  volumes  of  the  headspace  and  water  sample,  and  the  Henry's  law 
constant  of  the  particular  gas. 

30.  Use  of  the  headspace  analysis  technique  with  a  GC  is  much  more  com¬ 
plex  than  the  titrametric  or  colorimetric  techniques.  Experienced  laboratory 
personnel  are  required  to  operate  the  GC  and  prepare  the  samples.  However, 
for  hydrocarbon  tracer  gases,  this  technique  proves  to  be  one  of  the  best. 
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31.  Any  adequate  analysis  requires  that  the  accuracy  and  precision  of 
the  data  be  assessed.  Uncertainty  analysis  techniques  provide  the  tools  and 
methodology  to  perform  this  evaluation.  This  technique  can  be  used  to  evalu¬ 
ate  existing  data,  or  from  a  planning  perspective,  it  can  provide  guidance  on 
data  requirements  to  assure  that  high-quality  data  are  collected,  which  will 
provide  a  proper  foundation  for  drawing  conclusions.  For  in  situ  or  discrete 
sampling  methodologies,  the  uncertainty  in  the  collected  data  must  be  evalu¬ 
ated  for  confidence  to  exist  in  any  subsequent  conclusions.  Details  of  the 
analysis  technique  are  presented  in  Appendix  A  for  a  typical  data  collection 
program. 

32.  If  a  first-order,  second-moment  uncertainty  analysis  (Rindels  and 
Gulliver  1989)  is  performed  upon  the  measurements  that  go  into  determining  the 
transfer  efficiency,  the  uncertainty  in  the  efficiency  can  be  expressed  as: 


where 

Ue  -  total  uncertainty  in  E 

uncertainties  in  and  Cj  ,  respectively 

Be  -  bias  uncertainty  in  the  measurement  of  and  Cf 

Bc^  “  bias  uncertainty  in  C, 

Equation  10  shows  that  the  uncertainty  in  E  is  inversely  related  to  the  ini¬ 
tial  (upstream)  deficit.  The  typical  precision  uncertainty  for  well- 
calibrated  DO  meters  is  0.1  mg/i;  thus,  it  can  be  assumed  that  and  Wg^ 

are  0.1  mg/ Ji  at  the  95  percent  confidence  interval  (P  -  0.95).  The  variation 
in  time  of  DO  concentration  at  a  given  location  can  result  in  a  Wg^  or  Wg^ 
value  that  is  somewhat  larger.  The  last  two  terms  of  Equation  10  represent 


the  bias  uncertainties  that  need  to  be  considered.  These  uncertainties  repre¬ 
sent  potential  error  that  would  be  identical  in  every  measurement  at  that 
structure,  or  that  day,  in  terms  of  DO  meter  calibration  and  the  uncertainty 
in  the  value  of  saturation  concentration.  A  typical  bias  uncertainty  for  DO 
probe  calibration  would  be  ±0.1  (P  -  0.95).  In  experiments  on  river  water 
sampled  during  midwinter  in  Minnesota,  Rindels  (1990)  found  that  the  DO  satu¬ 
ration  concentration  was  98  percent  of  that  for  distilled  water,  with  an  un¬ 
certainty  of  ±2  percent.  Thus,  Rindels  found  that  Bg  -  0.02C,  .  During 
summer,  because  of  the  decreased  water  quality,  Bg  may  be  higher,  say  Bg  - 
O.OlCj  .  This  is  the  primary  source  of  uncertainty  in  many  measurements  of 
transfer  efficiency. 

33.  As  mentioned  at  the  beginning  of  this  section,  uncertainty  analysis 
can  be  used  in  several  types  of  measurements.  Consider  the  following  ques¬ 
tion,  which  would  be  posed  prior  to  a  field  study:  How  much  DO  deficit  must 
exist  to  accurately  assess  the  oxygen  transfer  characteristics  of  a  hydraulic 
structure?  Uncertainty  analysis  provides  the  framework  for  establishing  those 
minimum  criteria  for  in  situ  measurements  of  DO.  For  example,  if  the  oxygen 
saturation  value  Cg  is  8.0  mg/i,  the  expected  efficiency  is  50  percent  (Dj 

=  0.5Di)  ,  a  3  percent  bias  uncertainty  in  C,  exists,  and  it  is  decided  that 
the  uncertainty  in  E  must  be  less  than  10  percent  (U^  <  0.30E)  to  use  in 
this  analysis,  then  the  results  of  applying  Equation  10  are  that  D^  must  be 
greater  than  3.4  rag/f.  If  these  conditions  exist  for  this  study,  e.g.,  C,  - 
8.0  mg/i  ,  measured  E  =  0.50  ,  D^  -  3.4  mg/i  ,  then  the  measured  efficiency 
would  be  0.50  (50  percent)  ±0.05  (or  10  percent  of  E)  for  a  confidence  inter¬ 
val  of  95  percent.  If  the  upstream  deficit  is  less  than  2.5  mg/i,  then  the 
assessment  of  the  gas  transfer  characteristics  of  the  structure  should  not  be 
based  on  the  measurement  of  oxygen  concentrations.  An  alternative  method 
using  a  tracer  gas  should  be  considered. 

34.  After  a  field  study,  the  uncertainty  surrounding  the  observed  data 
should  be  evaluated,  as  detailed  in  Appendix  A,  so  that  some  measure  of  con¬ 
fidence  can  be  defined  about  the  study.  For  an  upstream  oxygen  concentration 
of  3  mg/i,  a  downstream  concentration  of  5.5  mg/i,  and  saturation  of  8  mg/i, 
the  efficiency  determined  with  these  data  from  Equation  10  is  given  as  E  -  0.5 
±  0.034  (P  =  0.95).  The  uncertainty  here  is  approximately  7  percent  of  the 
measured  efficiency.  Obviously,  the  uncertainty  in  E  becomes  smaller  if  the 
upstream  deficit  is  larger.  Hence,  it  is  easy  to  conclude  that  there  will  be 
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less  uncertainty  if  measurements  are  made  when  the  magnitude  of  the  upstream 
deficit  is  maximum. 

35.  The  measurement  uncertainty  technique  can  be  applied  to  any  type  of 
measurement,  if  the  uncertainty  of  each  part  can  be  defined.  In  situ  measure¬ 
ments  of  DO  to  determine  efficiency  were  used  as  an  Illustration,  but  all  the 
measurement  techniques  discussed  may  be  analyzed  in  this  fashion,  using 
Equations  Al,  A2 ,  and  A3. 
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PART  IV:  HYDRAULICS  AT  VARIOUS  STRUCTURES 


36.  Low-head  structures  range  widely  in  purpose,  size,  and  configura¬ 
tion.  Consequently,  the  hydraulics  that  are  encountered  at  the  variety  of 
structures  are  extremely  varied.  However,  the  processes  identified  in  Part  II 
are  still  the  major  influences  on  gas  transfer.  An  examination  of  the  hydrau¬ 
lic  conditions  at  "generic"  structures  can  reveal  the  impact  of  these  proc¬ 
esses  and  permit  qualitative  conclusions  to  be  drawn  regarding  the  relative 
gas  transfer  characteristics  of  these  structures. 

37.  A  committee  evaluation  of  gas  transfer  at  hydraulic  structures 
(American  Society  of  Civil  Engineers  (ASCE)  1991)  categorized  open-channel 
structures  into  four  groups:  (a)  free-surface  flows,  such  as  flow  in  a  chan¬ 
nel  or  on  a  spillway  or  ogee  crest  without  a  tailwater  plunge  pool  (Figure  3); 

(b)  sub-merged  flows,  such  as  discharge  under  a  submerged  gate  (Figure  4); 

(c)  free  jets,  such  as  flow  over  a  sharp-crested  weir  (Figure  5);  and 

(d)  transitional  flows,  where  free-surface  flows  or  jets  interact  with  a  pool 
of  water  resulting  in  plunging  flow  or  a  hydraulic  jump  (Figures  6  and  7), 

The  hydraulics  of  each  group  differ  greatly  resulting  in  significantly 
different  gas  transfer  characteristics.  An  understanding  of  the  hydraulics 
and  reaeratlun  character  of  these  groups  of  structures  can  permit 
extrapolation  to  more  complex  flow  situations. 

38 .  Free-surface  flows ,  such  as  shown  in  Figure  3 ,  generate  a  boundary 
layer  along  the  spillway  surface.  If  the  length  of  flow  is  sufficient,  then 
the  turbulence  of  this  boundary  layer  will  intersect  the  free  surface.  Under 
these  conditions,  air  is  captured  by  the  highly  turbulent  surface  and  sheared 
down  to  small  bubbles,  giving  the  flow  a  white  appearance.  As  the  entrained 
air  is  carried  along  the  spillway,  the  bubbles  are  transported  downward 
through  the  water  column  by  turbulence.  As  a  consequence  of  this  entrained 
air,  the  surface  area  available  for  gas  transfer  greatly  increases  compared  to 
spillway  flows  that  do  not  experience  free-surface  aeration.  Gulliver,  Thene, 
and  Rindels  (1990)  estimated  that  the  surface  area  of  the  bubbles  in  a  1-ft- 
deep  spillway  flow  can  provide  500  times  the  area  of  air/water  interface  pro¬ 
vided  by  the  free  surface.  At  overflow  crests  where  air  entrainment  occurred 
in  the  free-surface  flow,  Rindels  and  Gulliver  (1989)  observed  oxygen  transfer 
efficiencies  as  high  as  30  percent,  which  represents  approximately  half  of  the 
oxygen  absorption  occurring  at  the  project. 


23 


(after  ASCE  1991) 


(after  ASCE  1991) 
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Figure  6.  Transitional  flows  through  a  gated  spillway 

(after  ASCE  1991) 


Poulbto  oir  •nlroinfiMnl  ot 
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39.  Submerged  discharges  (Figure  4)  usually  reaerate  flow  throu^  oxy¬ 
gen  absorption  at  the  water  surface  of  the  turbulent  "boil"  in  the  stilling 
basin  of  the  structure.  Generally,  there  is  little  air  entrainment  as  long  as 
the  discharge  remains  submerged.  As  a  consequence,  the  reaeration  that  occurs 
in  these  hydraulic  conditions  is  relatively  small  with  oxygen  absorption  effi¬ 
ciencies  usually  below  10  percent  (Wilhelms  1988;  Thene,  Daniil,  and  Stefan 
1989).  However,  for  many  structures  like  the  one  shown  in  Figure  4,  with 
higher  discharges  there  is  sufficient  momentvim  in  the  jet  Issuing  from  under 
the  gate  to  push  the  tailwater  downstream  and  expose  the  jet  to  the  atmo¬ 
sphere.  This  situation  significantly  alters  the  hydraulic  characteristics  of 
the  structure.  The  flow  conditions  may  more  closely  resemble  those  shown  In 
Figure  6,  where  a  high-velocity  jet  interacts  with  a  pool  of  water  causing  a 
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hydraulic  Jump.  When  this  occurs,  large  volumes  of  air  are  usually  entrained 
at  the  velocity  discontinuity  between  the  jet  and  the  static  pool.  With  the 
entrained  air  bubbles  and  possibly  hydrostatic  pressure  on  those  bubbles, 
oxygen  absorption  significantly  Increases  to  efficiencies  as  high  as 
40  percent  (Wilhelms  1988;  Thene,  Daniil,  and  Stefan  1989). 

40.  Free  jets,  such  as  shown  in  Figure  5,  appear  to  be  similar  to  the 
free-surface  flow  on  a  spillway.  These  flows,  hoyever,  do  not  demonstrate  the 
free-surface  air  entrainment  that  results  from  the  full  development  of  the 
boundary  layer  at  the  spillway  surface.  Typically,  reaeratlon  is  accomplished 
at  this  type  of  structure  during  the  breakup  of  the  jet  when  it  collides  with 
*-.he  bottom  downstream.  If  the  free  jet  plunges  into  a  receiving  pool,  air 
entrainment  and  turbulent  mixing  contribute  to  increasing  gas  transfer. 
Further,  the  depth  of  the  plunge  pool  can  enhance  the  absorption  because  of 
the  increased  hydrostatic  pressure  on  the  entrained  air  bubbles.  Avery  and 
Novak's  (1978)  experiments  indicate  that  the  transfer  efficiency  increases 
with  a  tailwater  depth  up  to  0.6  times  the  drop  height.  Oxygen  absorption 
efficiencies  vary  widely,  but  for  low-head  overflow  weirs,  efficiencies  of  up 
to  70  percent  have  been  measured. 

41.  The  most  complex  situations  occur  in  open  channel  transitional 
flows,  which  are  typical  of  many  hydraulic  structures  with  stilling  basins. 
Much  of  the  time,  these  situations  are  combinations  of  several  of  the  hydrau¬ 
lic  conditions  described  in  the  preceding  paragraphs.  For  example,  this  may 
include  air-entrained  flow  over  a  crest  that  plunges  into  a  tailwater;  the 
hydraulic  conditions  when  discharges  from  a  gated  sill  sweep  the  tailwater 
downstream  and  expose  the  flow  to  the  atmosphere,  resulting  in  air  entrain¬ 
ment;  or  the  variety  of  flow  conditions  encountered  in  a  gated  conduit  outlet. 
Obviously,  the  hydraulic  action  of  these  flow  conditions  is  complex  and 
extremely  varied.  Consequently,  the  reaeration  characteristics  of  these  flows 
are  extremely  varied  with  oxygen  transfer  efficiencies  ranging  up  to  nearly 
100  percent. 
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PART  V;  FIELD  DATA  AND  PREDICTION  OF  REAERATION 


42.  Qualitative  comparisons  can  be  made  of  the  reaeration  that  occurs 
under  various  flow  conditions  by  considering  the  significance  of  the  major 
processes  that  affect  gas  transfer.  Air  entrainment  usually  results  in 
greater  gas  transfer,  as  does  plunging  aerated  flow.  Generally,  structures 
with  greater  differences  between  pool  and  tailwater  elevations  demonstrate 
greater  gas  transfer,  i.e.,  there  is  more  energy  available  to  cause  gas  trans¬ 
fer.  Through  observations  of  oxygen  absorption  or  tracer  gas  transfer  at 
these  structures,  the  reaeration  characteristics  of  the  structure  may  some¬ 
times  be  parametrically  described.  This  type  of  analysis  results  in  mathemat¬ 
ical  relationships  between  gas  transfer  and  the  measurable  parameters  that  are 
descriptive  of  the  significant  processes  at  work.  The  following  sections 
review  most  of  the  efforts  in  this  area,  including  data  presentation,  uncer¬ 
tainty  analysis,  and  an  evaluation  of  predictive  capability  of  various 
mathematical  descriptions . 


Field  Data 


43.  Parameter  measurement  and  subsequent  calculations  with  observed 
data  constitute  crucial  initial  steps  in  practically  all  engineering  applica¬ 
tions,  and  reaeration  studies  are  no  exception.  While  the  amount  of  informa¬ 
tion  needed  to  complete  a  thorough  gas  transfer  study  at  a  hydraulic  structure 
is  not  large,  measures  of  confidence  are  important  and  the  collection  of  data 
is  not  trivial  (Rindels  and  Gulliver  1989) .  Prediction  equations  also  require 
special  consideration  since  each  uses  a  unique  set  of  variables  and  special 
measurements  are  sometimes  necessary.  This  section  describes  how  data  used  to 
evaluate  the  predictive  equations  for  gas  transfer  at  hydraulic  structures 
have  been  gathered  and  analyzed. 

44.  The  entire  raw  data  set  used  in  the  analyses  reported  herein  is 
presented  in  Appendix  B,  which  was  created  by  reviewing  current  literature  and 
Includes  nvunerous  parameters,  all  in  SI  units.  The  following  parameters  are 
cited- 

a.  The  location  and  type  of  the  hydraulic  structure. 

b.  The  date  the  sampling  was  made. 

c.  The  flow  per  unit  width  over  the  structure. 
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The  difference  In  elevation  upstream  and  downstream  of  the 
structure  (the  head  loss) 

£.  The  depth  of  the  tallwater. 

The  gate  openings  If  applicable. 

g.  The  concentration  of  dissolved  oxygen  both  upstream  and 
downstream  of  the  structure. 

tl.  The  temperature  of  the  water. 

j,.  The  saturation  concentration  of  oxygen  In  the  water. 

J,.  The  barometric  pressure  at  or  near  the  site. 

45.  When  a  report  failed  to  Include  the  saturation  concentration  of 
oxygen  C,  the  value  was  found  using  the  equation  of  Hua  (1990): 


5355  ♦0.0226297*T+ 

0.01166  -  f. 

^llcc 

1  (12) 

1  T  J 

L  1 

—Tjr— JJ  C 

1 

where  T  Is  temperature  In  degrees  kelvln  and  Ccj.  Is  concentration  of 
chloride  Ions  In  grams  per  litre.  Chloride  concentration  is  normally  not  part 
of  the  data  collected  at  hydraulic  structures  because  for  concentrations  under 
1  g/i,  the  presence  of  chloride  has  little  Impact  on  saturation  concentra¬ 
tions.  However,  most  natural  streams  contain  some  chloride  Ions,  and  particu¬ 
larly  near  maritime  environments,  atmospheric  deposition  of  chloride  can  cause 
error  In  estimations  of  oxygen  saturation.  For  the  purposes  of  this  report, 
Cci  was  assumed  In  all  cases  to  be  1  g/i. 

46.  The  equation  of  Hua  expresses  the  saturation  concentration  of  oxy¬ 
gen  In  distilled  water  at  an  atmospheric  pressure  of  760  mm  mercury.  The 
value  was  adjusted  for  the  following: 

a.  Water  quality  by  estimating  the  value  of  C,  to  be  0.97  of  the 
distilled  water  value  (Rlndels  1990) . 

The  barometric  pressure  with 


C.  (P.un)  -  C.  (760 


mm  Hg 

®  '760 


(13) 


where  Is  the  reported  pressure  in  millimetres  of  mercury.  If 

barometric  pressure  was  unavailable,  then  the  saturation  concentration  was 
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adjusted  for  the  elevation  of  the  structure  above  sea  level  by  employing  the 
following  equation: 


Elevation  Correction 


ri 


(0.006SN) 

- 


TS.25 


(14) 


where  M  Is  elevation  above  mean  sea  level  In  metres. 

47.  The  processed  data  for  comparison  with  the  predictive  equations  are 
presented  In  Appendix  C.  This  appendix  lists  several  of  the  previous  param¬ 
eters  that  are  used  In  predictive  eqxiations.  In  addition  to  values  of  transfer 
efficiency  (Equation  8)  and  uncertainty  (Equation  10) .  The  efficiencies  pro¬ 
vided  by  the  literature  were  commonly  calculated  at  a  variety  of  different 
temperatures.  For  comparison,  these  efficiencies  were  Indexed  to  a  common 
temperature  (20  *C)  using  the  equation  of  Gulliver,  Thene,  and  Rlndels  (1990): 

where  E20  and  Ex  are  gas  transfer  efficiencies  at  20  and  at  T  *’C  and  fx 

is  given  by  the  equation: 

fx  -  1.0  +  0.02103(t  -  20)  ♦  8.261  *  10'5(t  -  20f 

where  T  Is  temperature  In  degrees  Celsius. 

48.  The  uncertainty  listed  In  Appendix  C  Is  the  cumulative  uncertainty 
to  the  95  percent  confidence  Interval  associated  with  the  gas  transfer  effi¬ 
ciency.  It  Is  calculated  by  Equation  10  or  a  similar  equation,  using  the 
principles  outlined  In  Appendix  A.  The  bias  uncertainty  In  saturation  concen¬ 
tration  accounts  for  the  fact  that  the  saturation  for  DO  Is  not  precisely 
established  In  field  situations.  Rlndels  (1990)  found  that  during  midwinter 
this  bias  was  approximately  ±2  percent.  For  the  general  data  sets,  however, 
and  accounting  for  the  degraded  water  quality  of  the  growing  season,  the  bias 
In  C,  was  taken  to  be  ±3  percent.  The  calibration  bias  for  DO  meters  was 
taken  to  be  ±1  percent. 
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49.  The  maxlmuni  acceptable  value  of  uncertainty  U  was  set  at  0.25  for 
values  of  transfer  efficiency  to  be  included  in  the  final  data  set  for  compar¬ 
ison  with  the  predictive  equations.  Most  of  the  literature  did  not  contain 
uncertainty  analyses  of  the  field  data;  therefore,  appropriate  approximations 
were  made  in  calculating  U  . 


Predictive  Equations 

50.  As  discussed  at  the  beginning  of  this  section,  predictive  equations 
that  have  been  developed  by  various  researchers  for  gas  transfer  at  hydraulic 
stiructures  generally  use  physical  parameters  of  the  structure  or  flow  condi¬ 
tions,  l.e. ,  Froude  number,  depth  of  tailwater,  discharge,  etc.,  to  estimate 
reaeratlon  efficiency.  These  equations  are  all  more  or  less  empirical,  with 
intuition  being  used  to  specify  the  independent  variables,  and  regressions 
used  to  determine  the  constants  in  the  equation.  Each  equation  yields  a 
unique  oxygen  transfer  efficiency,  and  it  is  difficult  to  know  the  correct 
efficiency  for  a  different  site  without  performing  a  comprehensive  field  mea¬ 
surement  investigation.  The  observed  data  discussed  in  the  previous  section 
and  presented  in  Appendix  B  are  used  in  a  comparison  with  the  predictive  equa¬ 
tions.  Predicted  and  measured  oxygen  transfer  efficiencies  are  compared  to 
determine  the  performance  of  the  various  predictive  equations.  It  should  be 
recognized  that  predictive  equations  are  reliable  only  as  first  approxima¬ 
tions,  and  that  they  should  be  used  carefully  only  when  more  accurate  and 
detailed  studies  are  not  feasible. 

51.  The  equations  investigated  and  the  references  in  which  they  may  be 
found  are  listed  in  Table  1.  A  more  thorough  explanation  on  the  development 
of  these  equations  is  presented  in  Appendix  D. 

52.  Each  of  these  equations  is  derived  from  a  different  data  set,  with 
coefficients  adjusted  to  fit  the  specific  data  of  the  researcher.  Because  the 
size  of  most  data  sets  is  relatively  small,  large  deviations  between  predicted 
values  of  gas  transfer  efficiency  and  measured  values  at  different  sites  are 
common.  Appendix  D  shows  plots  of  predicted  versus  measured  oxygen  transfer 
efficiency,  which  illustrate  this. 

53.  Each  of  the  equations  in  Table  1  was  used  to  predict  an  oxygen 
transfer  efficiency  for  all  of  the  acceptable  measured  data  (U  <  0.25)  found 
in  Appendix  C.  This  was  performed  separately  for  each  of  the  four  categories 
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Table  1 

Predictive  Equations  of  Gas  Transfer  at  Hydraulic  Structures 


Hydraulic 

Structure 


Sharp- 
Crested 
Heirs 
(Figure  5) 


Predictive  Equation 


fox  ^  *  l.SHe)  >  1.2  (n)  and  q  i  235 


Source 


Avery  and 
Novak 
(1978) 


Nakasone 

(1987) 


E,o  •  1  -  axp -[0.0861  ^  ♦  1.5He)*-*  “q*-***H*-»^®] 


for  (h  *  l.SH,)  1.2  (m)  and  q  >  235 


Bjo  ■  1  -  exp -[5. 39  (h  ♦  1.5Ho)»-”q-«  >"H®  "•) 


fox  (h  *  1.5He)>l.2  (m)  and  q  >  235 


B„  »  1  -  exp-[5.92  (h  ♦  l.5H„)®-“*q-®  “>H®  ”•] 


fox  ^  *  1.5H,)  s  1.2  (b)  and  q  i  235 


B,,  -  1  -  exp -[0.0785  (h  ♦  1.5H„)^-“q»-*”H»  ”») 


xp  -(0.156  »*• 


(Continued) 


Thene 

(1988) 
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Table  1  (Continued) 


Hydraulic 

Structure 


Ogee  Crest 
Spillways 
(Figure  3) 


Gated  Sills 
and  General 
Hydraulic 
Structures 
(Figure  4) 


Gated 
Conduits 
(Figure  7) 


Predictive  Equation 


1  +  1.005  X  10-*  rf  **  R*  “  fl  -  0.6  exp 


(->•’  Dll' 


B,.  •  1  -  exp  f — -0.2625h  _  0.2034H 

^  I  1  ♦  0.2153q 


E,4  -  1  -[exp(-0.5249h)]®  *"* 


E,o  "  1  -  exp  (-0.1772h) 


Bj,  «  exppo. 00857884.^  -  0.188 


E,,  -  1  -  exp  (-0.1476b) 


Source 


Thene 

(1988) 


Holler 

(1970) 


Rlndels  and 
Gulliver 
(1991) 


Foree 

(1976) 


Tsivoglou 
and  Wallace 
(1972) 


Freul  and 
Holler 
(1969) 


Wilhelms 

(1988) 


Wilhelms 
and  Smith 
(1981) 


where 


Fj  -  Froude  number  of  the  jet 


(Continued) 
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Table  1  (Concluded) 


Hydraulic 

Structure 

Predictive  Equation 

Source 

P  ,  (2g)°-2ShO  » 

J  TTT3 


g  -  acceleration  due  to  gravity,  m/sec^ 
h  »  head  loss,  m 
q  -  specific  discharge,  m^/sec/m 
R  -  Reynolds  number 


R  = 


q 


1/  -  kinematic  viscosity 
He  -  critical  water  depth  on  the  weir,  m 
H  -  tailwater  depth,  m 

Nf  -  Froude  number  at  the  point  of  impact 


Nf 


S 


(2gh)°” 


S  ~  gate  submergence,  m 


of  hydraulic  structures  (an  ogee  crest  spillway,  gated  sill,  weir,  or  gated 
conduit)  regardless  of  which  type  of  structure  the  predictive  equation  was 
developed  for.  The  statistical  results  are  given  in  Table  2.  The  standard 
error  compares  measured  and  predicted  oxygen  transfer  efficiencies  by  the 
equation 


SE 


(E^  -  Ep)^ 

n 


(17) 


where  and  Ep  are  measured  and  predicted  oxygen  transfer  efficiency, 

respectively.  The  standard  error  is  approximately  equal  to  the  uncertainty  U 
to  the  68  percent  confidence  interval.  A  68  percent  confidence  interval  may 
be  interpreted  as  meaning  that  two  out  of  three  predicted  values  will  be 
within  the  boundaries  of  ±U.  An  uncertainty  to  the  95  percent  confidence 
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Table  2 

Summary  of  Statistical  Results 

Formula  Title 

Standard  Error 

Ogee 

Gated 

Sill 

Weir 

Gated 

Conduit 

Avery  and  Novak  (1978) 

0.282 

0.458 

0.166 

0.340 

Thene's  Adjustment  to  Avery  and 
Novak  (1988) 

0.297 

0.451 

0.170 

— 

Preul  and  Holler  (1969) 

0.647 

0.141 

0.615 

0.690 

Thene  (1988) 

0.302 

0.324 

0.174 

0.420 

Nakasone  (1987) 

0.267 

0.487 

0.172 

— 

Tsivoglou  and  Wallace  (1972) 

0.290 

0.406 

0.183 

0.320 

Foree  (1976) 

0.285 

0.612 

0.271 

0.358 

Rindels  and  Gulliver  (1991) 

0.160 

0.463 

0.210 

— 

Holler  (1970) 

0.327 

0.296 

0.205 

0.339 

Wilhelms  (1988) 

0.227 

0.247 

0.360 

— 

Wilhelms  and  Smith  (1981) 

0.322 

0.355 

0.212 

0.312 

Interval  (such  as  that  used  for  the  measurements)  may  be  approximated  as 
double  the  uncertainty  of  the  68  percent  confidence  interval,  or: 


U 


E(P.0.95) 


2U, 


E(P.0.68) 


(18) 


where  P  is  the  probability  of  the  confidence  interval  on  efficiency. 

54.  Using  the  method  described  herein,  the  prediction  of  a  transfer 
efficiency  will  result  in: 


E 


E  ±  U, 


E(p 


0.68) 


(19) 


where 


E  -  true  value  of  transfer  efficiency  that  the  given  equation  is 
designed  to  predict 
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E  ••  the  best  estimate  of  the  true  value  of  transfer  efficiency 
that  is  predicted  by  the  equation 

^B(p>4.6S)  -  an  uncertainty  to  the  68  percent  confidence  interval 
(p  -  0.68)  for  the  prediction 

55.  While  specific  equations  may  provide  the  best  approximation  for 
reaeratlon  efficiencies  at  individual  sites,  the  equation  that  minimizes 
predicted  uncertainty  provides  the  best  estimate  for  gas  transfer  at  a  generic 
hydraulic  structure.  The  results  show  that  the  equation  of  Avery  and  Novak 
(1978)  best  estimates  the  transfer  efficiency  for  weirs.  Rindels  and 
Gulliver's  (1989)  equation  is  most  accurate  for  gated  and  tingated  ogee  spill¬ 
ways.  The  equation  of  Preul  and  Holler  (1969)  predicts  reaeration  at  gated 
sills  most  accurately.  The  equation  of  Wilhelms  and  Smith  (1981)  best 
describes  the  reaeration  of  gated  conduits,  although  none  of  the  equations 
predicted  the  process  accurately.  This  is  probably  because  of  the  tranr- 
sltlonal  characteristics  of  the  flow  in  the  conduit  (high-velocity  open  chan¬ 
nel  flow  to  plunging  flow  or  hydraulic  Jump  open  channel  versus  closed  conduit 
flow),  and  the  lack  of  consistency  between  structures. 
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PART  VI:  SUMMARY  AND  CONCLUSIONS 


56.  Gas  transfer  or  reaeratlon  is  considered  a  first-order  reaction 
process  and  can  be  conveniently  described  by 


The  Important  physical  processes  that  Impact  gas  transfer  are  Included  In  this 
formulation  in  the  following  manner:  The  effect  of  turbulent  mixing  Is  re¬ 
flected  by  the  liquid  film  coefficient  ki,  .  The  Impact  of  surface  area 
available  for  gas  transfer,  which  must  Include  the  surface  area  of  entrained 
air  bubbles,  Is  represented  by  the  specific  area  term  A/V  .  Enhanced  gas 
absorption  due  to  the  effects  of  hydrostatic  pressure  on  plunging  aerated  flow 
Is  Included  as  a  pressure  modification  of  the  saturation  concentration  C,  . 
The  time  of  contact  over  which  gas  transfer  can  occur  Is  t  . 

57,  The  unknowns  in  this  equation  often  dictate  that  reaeratlon  mea¬ 
surements  be  conducted  to  determine  the  gas  transfer  efficiency.  The  transfer 
efficiency  can  then  be  related  to  the  physical  processes  through  empirical 
relationships  and  regression  analyses.  Several  alternatives  are  available  for 
measurement  of  gas  transfer.  Direct  In  situ  measurement  of  dissolved  oxygen 
with  polarographlc  probes  Is  usually  the  most  convenient.  However,  based  on 
uncertainty  analysis,  a  DO  deficit  of  at  least  2.5  mg/i  (when  saturation  Is 
8.0  mg/i)  is  required  for  accurate  analysis.  Often,  the  DO  is  not  sufficient¬ 
ly  low  to  permit  an  accurate  analysis.  When  this  Is  the  situation,  alterna¬ 
tive  gas  transfer  tracers  must  be  used  in  measurements  for  gas  transfer.  In 
situ  methane  gas  shows  the  highest  potential  for  general  application  In 
reaeratlon  field  studies. 

58,  Generally,  low-head  structures  can  be  categorized  into  four  groups: 
(a)  free-surface  flows,  such  as  flow  in  a  channel  or  on  a  spillway  or  ogee 
crest  without  a  tallwater  plunge  pool;  (b)  submerged  flows,  such  as  discharge 
under  a  submerged  gate;  (c)  free  Jets,  such  as  flow  over  a  sharp-crested  weir; 
and  (d)  transitional  flows,  where  free-surface  flows  or  jets  Interact  with  a 
pool  of  water  resulting  in  plunging  flow  or  a  hydraulic  jump.  The  hydraulics 
of  each  group  differ  dramatically  and,  consequently,  the  gas  transfer  charac¬ 
teristics  are  significantly  different.  However,  an  understanding  of  the 
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hydraulics  and  the  resulting  gas  transfer  characteristics  of  each  type  of 
structure  can  permit  limited  extrapolation  to  other  hydraulic  structures.  The 
observations  reported  herein  show  that  gated  and  ungated  ogee  crests  demon¬ 
strate  gas  transfer  efficiencies  up  to  nearly  100  percent.  Submerged  flows 
without  air  entrainment  can  result  In  efficiencies  up  to  40  percent.  Sharp- 
crested  overflow  weirs  or  jets  demonstrate  efficiencies  up  to  70  percent. 
Transitional  flow  conditions,  which  Include  several  or  all  of  the  other  types 
of  flow  conditions,  have  shown  reaeratlon  efficiencies  of  up  to  nearly 
100  percent. 

59.  Several  sources  of  observed  data  were  assembled  and,  after  screen¬ 
ing  through  an  uncertainty  analysis,  are  Included  In  the  appendices  of  this 
report.  Predictions  from  eleven  equations  that  describe  gas  transfer  effi¬ 
ciency  at  various  types  of  structures  were  compared  to  this  comprehensive 
database.  Although  there  are  inadequacies  In  the  models  currently  used  for 
prediction  of  gas  transfer,  the  equations  listed  In  Table  3  were  the  most  suc¬ 
cessful  for  given  types  of  structures. 


Structure 

Type 


Ogee 

Crests 


Gated 

Sills 


Sharp- 

crested 

Weirs 


Gated 

Conduit 

Outlets 


Table  3 

Suggested  Predictive  Equations 


Predictive  Equation 


-0.262Sh -  0.2034H 

1  ♦  0.2153q 


Ejo  -  1  - 


(l  ♦  666N**  ”) 


Standard 

Error 

Source 

0.16 

Rlndels 

and 

Gulliver 

(1991) 

0.14 

Pruel  and 
Holler 
(1969) 

0.17 

Avery  and 
Novak 
(1978) 

0.31 

Wilhelms 
and  Smith 
(1981) 
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60.  Standard  error  estimations  for  the  equations  of  Preul  and  Holler 
(1969)  and  Rlndels  and  Gulliver  (1991)  may  be  too  low.  Much  of  the  sill  and 
ogee  spillway  database  Included  In  this  report  was  used  In  the  regression  of 
the  two  equations,  and  the  paucity  of  accurate  alternative  data  sources  close¬ 
ly  tailors  these  equations  to  this  evaluation.  While  the  compilation  of  accu¬ 
rate  data  Is  not  simple,  further  field  Investigations  are  needed  to  define  the 
performance  of  predictive  equations. 

61.  In  their  efforts  to  predict  the  Impact  of  hydraulic  structures  on 
levels  of  dissolved  gases  In  river  systems,  many  researchers  sought  to  create 
equations  that  provide  accurate  estimates  of  transfer  efficiency.  Unfortu¬ 
nately,  the  development  of  each  equation  was  limited  by  the  size  of  the  data¬ 
base  that  a  researcher  had  available  and  the  difficulty  encountered  in  deriv¬ 
ing  a  theory  from  first  principles.  As  a  result,  many  equations  are  useful 
only  for  specific  types  of  structures  under  particular  conditions,  l.e.,  the 
domain  of  accuracy  for  the  parameters  used  In  the  equation  Is  relatively 
small.  The  relatively  large  deviations  between  the  measured  transfer  effi¬ 
ciencies  and  the  values  computed  from  the  predictive  equations  are  understand¬ 
able  because  of  the  size  and  diversity  of  the  observed  data.  Consequently,  a 
large  uncertainty  Is  associated  with  the  prediction  of  an  oxygen  transfer 
efficiency  from  these  equations.  Field  measurements  are  still  the  most  con¬ 
sistent  means  of  determining  oxygen  transfer  characteristics. 
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APPENDIX  A:  TYPICAL  MEASUREMENT  UNCERTAINTY  ANALYSIS  FOLLOWING 
THE  EVALUATION  TECHNIQUE  OF  RINDELS  AND  GULLIVER  (1989) 


A1 


1.  Measurement  uncertainty  inadvertently  occurs  due  to  the  inaccuracy 
of  the  instruments  or  of  the  measurement  or  because  of  operator  error.  In  any 
measurement  procedure,  all  significant  uncertainties  should  be  quantified,  to 
provide  a  means  of  evaluating  the  quality  of  the  procedure  and  resulting  data. 
During  the  course  of  the  measurements,  sources  of  uncertainty  should  be  iden¬ 
tified.  For  example,  in  the  chemical  titrametric  technique  for  determining 
dissolved  oxygen  (DO)  concentration,  the  purity  of  the  chemical  titrant  in  the 
Uinkler  technique  and  the  minimum  gradations  on  the  titration  buret  will 
contribute  to  the  uncertainty  of  the  measured  DO  concentrations  and  to  the 
uncertainty  in  transfer  efficiency.  In  the  use  of  polarographic  probes  for 
measuring  DO  concentration,  uncertainty  occurs  because  of  the  accuracy  and 
calibration  of  the  instrument. 

2.  From  Equations  8  and  9  in  the  main  text,  transfer  efficiency  E* 
was  defined  as 


E  =  1  - 


(Al) 


where  Cg  ,  Ci  ,  and  Cf  are  the  saturation  and  upstream  (initial),  and 
downstream  (final)  DO  concentrations,  respectively.  The  total  uncertainty  in 
the  transfer  efficiency  Ue  is  inherent  in  each  measurement.  By  definition, 
the  total  uncertainty  of  any  measurement  is  a  combination  of  precision  (ran¬ 
dom)  uncertainty  introduced  when  measurements  are  repeated,  and  bias  (syste¬ 
matic)  uncertainty,  or  possible  error  that  would  affect  each  measurement  in 
the  same  manner  (Abernethy,  Benedict,  and  Dowdell  1985).**  The  most  common 
technique  of  analyzing  measurement  uncertainties  is  a  first  order-second 
moment  analysis  (Kline  1985).  With  this  uncertainty  technique,  the  total 
uncertainty  of  transfer  efficiency  Ue  can  be  expressed  as 

U|  =  Wf  +  b|  (A2) 


*  For  convenience,  symbols  and  unusual  abbreviations  used  in  this  appendix 
are  listed  and  defined  in  the  Notation  (Appendix  E) . 

**  References  cited  in  this  appendix  can  be  found  in  the  References  at  the 
end  of  the  main  text  of  the  report. 
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where 


U£  -  precision  uncertainty  in  the  transfer  efficiency 
&£  -  bias  uncertainty  in  the  transfer  efficiency 

3.  The  precision  uncertainty  U£  is  a  conbination  of  the  sampling  un¬ 
certainties  that  result  from  determining  and  and  is  mathematically 

described  by  the  following  formula; 


I 


1(1  -  e)  W, 


Cl 


(A3) 


"e  and  are  the  uncertainties  associated  with  the  upstream  and 

.  stream  DO  measurements,  respectively.  According  to  Abernethy,  Benedict, 
and  Dowdell  (1985),  and  can  be  calculated  as  follows: 


t*  Or 


(A4) 


(A5) 


where 

t*  -  Student's  t-value  corresponding  to  n^  -  1  or  ng  -  1  degrees 
of  freedom  and  a  95  percent  confidence*^  Interval,  as  given  in 
Table  Al 

oq  -  standard  deviation  of  the  upstream  and  downstream  measurements, 
respectively 

n(;^,nc^  -  number  of  upstream  and  downstream  measurements,  respectively 

For  the  given  example  of  analysis  with  the  titrametrlc  method,  or 

should  be  assigned  a  minimum  value  of  0.05  mg/i,  which  is  the  uncertainty 
associated  with  reading  the  buret  used  for  titrations.  If  or  ,  as 

calculated  with  Equation  A3,  is  greater  Chan  0.05  mg/i,  then  their  respective 
calculated  values  would  be  used  in  determining  the  precision  uncertainty  W£ 
from  Equation  A3.  For  measurements  with  a  DO  meter,  or  Wc^  should 


Table  Al 

Student's  t-Values  at  a  95 

Percent  Confidence  Interval 

for  Various  Degrees  of  Freedom  (n^  -  1)  | 

Degree  of 

Student's 

Freedom 

t-Value 

1 

12.7 

.  2 

4.30 

3 

3.18 

4 

2.78 

5 

2.57 

6 

2.45 

7 

2.37 

8 

2.31 

9 

2.26 

10 

2.28 

15 

2.13 

20 

2.09 

30 

2.04 

have,  as  minimum  values,  the  stated  accuracy  of  the  DO  meter,  which  is  assumed 
to  be  the  uncertainty  at  the  95  percent  confidence  interval.  For  most  probes, 
this  is  usually  1  percent  of  full  scale  or  0.1  mg/i,  whichever  is  larger. 

4.  Two  bias  uncertainties  should  also  be  included  in  the  uncertainty: 
(a)  either  the  purity  of  the  chemical  titrant  in  the  Winkler  technique  or  the 
uncertainty  in  the  calibration  of  the  DO  meter  and  (b)  the  uncertainty  in  the 
value  of  saturation  concentration.  The  supplier  of  the  sodium  thiosulfate 
titrant  guarantees  the  purity  to  within  ±1  percent.  This  value  could  be 
assumed  to  be  the  uncertainty  at  the  95  percent  confidence  interval.  Proper 
air  calibration  of  a  DO  probe,  considering  air  pressure  and  relative  humidity, 
usually  has  a  typical  bias  uncertainty  of  ±0.1  (P  -  0.95).  The  saturation 
concentration  of  the  water  at  a  given  hydraulic  structure  was  estimated  using 
the  locally  measured  atmospheric  pressure,  before  and  after  sampling,  and  the 
water  temperature  measured  at  the  hydraulic  structure.  At  several  hydraulic 
structures,  Rlndels  (1990)  found  the  saturation  concentrations  to  be  less  than 
the  published  values  for  a  distilled  water  at  a  given  temperature.  He 
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estimated  the  actual  saturation  value  to  be  98  percent  of  the  distilled  water 
value  with  an  uncertainty  of  ±2  percent  for  river  water  collected  during 
winter.  Summertime  values  of  the  same  water  have  been  assumed  to  be 
97  percent  of  the  distilled  water  value  with  an  uncertainty  of  ±3  percent.  In 
situ  saturation  measurement  techniques,  currently  being  developed,  should 
reduce  this  uncertainty  considerably. 

5.  These  two  bias  uncertainties  are  added  to  give  the  total  bias  uncer¬ 
tainty  in  E  as  follows: 


*  (is  “'J  ■  [sr^J  * 


where  ,  the  bias  uncertainty  in  C,  ,  ranges  from  0.02  C,  to  0.03  C,  ; 

and  Bp  is  the  bias  uncertainty  in  DO  concentration  due  to  the  titrant  or  DO 
probe  calibration 


Be  ■  a  C, 


where  a  equals  0.02  for  the  titrant  and  0.01  for  the  probe. 

6.  Substituting  the  expressions  for  precision  uncertainties  (Equa¬ 
tions  A4  and  A5)  into  Equation  A3  gives  a  mathematical  description  of  the 
total  precision  uncertainty  In  E  .  Substituting  the  expressions  for  bias 
uncertainty  due  to  saturation  concentration  and  titrant  impurity  or  probe 
calibration  into  Equation  A6  gives  a  mathematical  description  of  the  total 
bias  uncertainty  in  E  .  Combining  these  in  Equation  A2,  the  uncertainty  in 
the  determination  of  transfer  efficiency  is  described  by 
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7 .  An  examination  of  Equations  A8  and  A9  clearly  shows  that  the 
accuracy  of  the  transfer  efficiency  is  strongly  dependent  on  the  difference 
between  saturation  and  the  measured  upstream  DO  concentration  (the  upstream 
oxygen  deficit).  Hence,  the  recommendation  was  made  in  paragraph  33  for 
minimum  upstream  deficits  for  acceptable  uncertainty  in  transfer  efficiency 
with  measurements  of  DO.  In  fact,  larger  oxygen  deficits  result  in  a  reduced 
measurement  uncertainty. 
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OBSERVED  DATA  TAKEN  FROM  LITERATURE 
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1.  For  some  references  all  the  information  that  was  needed  to  calculate 
Important  data  was  not  included  in  the  publications.  In  these  cases  allow¬ 
ances  had  to  be  made  to  extrapolate  the  information  from  different  sources. 

For  instance,  if  the  barometric  pressure  was  not  included,  then  (on  data  taken 
within  the  United  States)  the  nearest  National  Weather  Service  (NWS)  recording 
station  was  found  and  barometric  pressure  was  taken  from  NWS  publications. 

The  pressure  was  then  converted  to  the  elevation  of  the  structure,  because  NWS 
pressure  readings  have  all  been  converted  to  mean  sea  level  before  bei’' 
published. 

2.  Barometric  pressure  was  not  reported  for  these  structures  and  was 
gathered  from  NWS  records . 

Dresden  Island  Dam,  Joliet,  IL 

Brandon  Rock  Dam,  Joliet,  IL 

Starved  Rock  Dam,  Peoria,  IL 

Arkabutla  Dam,  Arkabutla,  MS 

Crooked  Creek  Dam,  Ford  City,  PA 

East  Branch  Dam,  Johnsonberg,  PA 

East  Lynn  Dam,  Wayne,  IN 

Enid  Dam,  Enid,  MS 

Grayson  Dam,  Grayson,  KY 

Grenada  Dam,  Grenada,  MS 

Mississinewa  Dam,  Peru,  IN 

Salamonie  Dam,  Largo,  IN 

Sardis  Dam,  Sardis,  MS 

Dams  listed  in  Butts  and  Evans  (1983)* 

3.  Although  pressure  was  not  reported  for  some  structures,  the  satura¬ 
tion  concentration  given  in  the  references  was  assumed  to  be  the  best  avail¬ 
able  estimate  for  the  following  structures: 

Borgharen  Weir,  Netherlands 
Lith  Weir,  Netherlands 
Cascade ,  Netherlands 

All  the  other  structures  in  Nakasone  (1979) 


*  References  cited  in  this  appendix  can  be  found  in  the  References  at  the  end 
of  the  main  text  of  the  report. 
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All  of  the  structures  listed  by  Foree  (1976) 


4.  An  "equivalent"  specific  discharge  for  all  of  the  gated  conduits  was 
calculated  from  hydraulic  radius  considerations  by  using  Fig.  6-5  of  Chow 
(1959)  for  a  fully  developed  velocity  profile  in  a  partially  full  circular 
conduit.  Then: 


q  - 


Q 


(Bl) 


where 

q  -  specific  discharge  through  the  structure 
Q  -  total  discharge  through  a  structure 
Rh  -  the  hydraulic  radius  of  the  flow  in  the  conduit. 

For  the  Japanese  dams  listed  by  Nakasone  (1979),  a  mean  elevation  of  550  ft 
was  assumed,  and  a  typical  air  pressure  at  mean  sea  level  was  determined  by 
inspecting  NWS  records.  The  measurement  uncertainty  associated  with  all  of 
the  structures  was  greater  then  Ug  -•  0.25,  so  none  were  included  in  the 
comparison  with  predicted  equations. 
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APPENDIX  D:  DESCRIPTION  OF  PREDICTIVE  EQUATIONS  AND  PLOTS  OF 
PREDICTED  VERSUS  MEASURED  OXYGEN  TRANSFER  EFFICIENCIES 


D1 


Avery  and  Novak  Equation 


1.  The  Reynolds  number  and  the  Froude  number  of  a  flow  are  used  by 
Avery  and  Novak  (1978)*  to  predict  the  deficit  ratio  specifically  for  weirs 
as  follows: 


ri5  -  1  -  0.64  X  10-*  Fr]  ”^  R°-533  (Dl) 

where  r^j  -  the  deficit  ratio  at  15  *C** 

Ftj  -  the  Froude  number  of  the  jet,  given  by  the  following 
equation; 


r-.  -  (2g)°  ” 

J  ^ - 

g  -  the  acceleration  due  to  gravity 
h  -  the  head  across  the  structure 
q  -  the  discharge  over  the  structure 
R  -  the  Reynolds  number  of  the  jet: 


1/  -  the  kinematic  viscosity 

The  deficit  ratio  r  can  then  be  related  to  the  efficiency  E  of  the 
structure  by: 


E 


(D2) 


*  References  cited  in  this  appendix  can  be  found  in  the  References  at  the 
end  of  the  main  text  of  the  report. 

**  For  convenience,  symbols  and  unusual  abbreviations  used  in  this  appendix 
are  listed  and  defined  in  the  Notation  (Appendix  E) . 


D3 


Also,  using  the  relation  of  Gulliver,  Thene,  and  Rlndels  (1990): 


1 

E20  -  1  -  (1  ' 


(D3) 


where  Ej  -  gas  transfer  efficiency  at  a  given  temperature  T 
fj  -  Indexing  coefficient  given  by: 

fi  -  1.0  +  0.02103  (T  -  20)  +  8.261  *  lO"*  (T  -  20)^ 

which  yields  an  efficiency  Indexed  to  20  *C  of: 


-20 


1  +  0 . 64  X  10' 


j^0,533 


1.1149 


(D4) 


The  results  of  the  Avery  and  Novak  equation  for  the  field  data  are  presented 
in  Figures  D1-D4. 


Adjustment  of  Thene  to  the  Averv  _and_ Novak  Eauatlon 


2,  Thene  (1988)  felt  that  the  accuracy  of  the  Avery  and  Novak  equation 
could  be  improved  by  considering  the  tallwater  depth  of  the  flow.  Absorption 
Is  expected  to  Increase  as  the  depth  and  contact  time  of  entrained  bubbles 
Increase.  Thene  created  a  function  that  varied  with  respect  to  the  depth  of 
the  tallwater  pool  H  ,  and  the  head  across  the  structure  h  .  The  adjustment 
of  Thene  to  Avery  and  Novak  takes  the  following  form; 


-15 


1  -  1.005  X  10-5 


p^2.08  £^0,63 


L  -  0.6  exp  |-3.7 


(05) 


which  can  then  be  converted  to  an  efficiency  at  20  °C  by  Equations  D2  and  D3 
to  yield: 


-20 


1 

'  1  +  1.005  X  10*5  FrJ  ®®  R°-55 

k 

L  -  0.6  exp  - 

k  k 

3.7  j 

^1 

1.1149 


(D6) 


D4 


Predicted  E 


•D  0.6 

(U 


0)  0.4 


0.0  0.2  0.4  0.6  0.8  1.0  1.2  1 

Measured  E 

Figure  Dl.  Avery  and  Novak  equation  for  ogee  spillways. 
Predicted  values  of  gas  transfer  versus  measured  values 
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Figure  D2.  Avery  and  Novak  equation  for  gated  conduits. 
Predicted  values  of  gas  transfer  versus  measured  values 


h 


Measured  E 

Figure  D3.  Avery  and  Novak  equation  for  weirs. 
Predicted  values  of  gas  transfer  versus  measured 

values 


Measured  E 

Figure  D4.  Avery  and  Novak  equation  for  gated 
sills.  Predicted  values  of  gas  transfer  versus 
measured  values 

D6 


Thene  reasoned  that  as  the  depth  of  the  tallwater  pool  H  increases ,  less 
flow  reaches  the  bottom  and  the  Importance  of  H  for  gas  transfer  diminishes. 
Because  of  this,  the  fall  height  of  the  structure  h  becomes  more  crucial  in 
estimating  the  amount  of  absorption  that  will  occur.  The  coefficients  of  the 
roots  in  the  original  Avery  and  Novak  equation  were  altered  to  improve  the 
statistical  accuracy  of  the  adjusted  formula.  Figures  D5-D7  show  Thene 's 
adjustment  to  the  Avery  and  Novak  equation. 


Thene  Equation 


3.  Thene  (1988)  extended  the  work  of  Elsawy  and  McKeogh  (1977)  into  a 
predictive  gas  transfer  equation  for  weirs.  For  a  rectangular  jet,  Ervine  and 
Elsawy  (1975)  found  the  air  entrainment  rate  to  predicted  by: 


(D7) 
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Q_ 
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Figure  D5.  Thene 's  adjustment  to  the  Avery  and  Novak  equation 
for  ogee  spillways.  Predicted  values  of  gas  transfer  versus 

measured  values 


D7 


Measured  E 

Figure  D6.  Thene's  adjustment  to  the  Avery  and 
Novak  equation  for  weirs.  Predicted  values  of 
gas  transfer  versus  measured  values 


0.0  0.2  0.4  0.6  0.8  1.0 
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Figure  D7.  Thene's  adjustment  to  the  Avery  and 
Novak  equation  for  gated  sills.  Predicted 
values  of  gas  transfer  versus  measured  values 


D8 


where 


Q,  and  Q*  -  air  and  water  discharges,  respectively 

b,  p,  t,  and  v  -  width,  perimeter,  thickness,  and  velocity  of  the  jet, 
respectively 

h  -  fall  height 

Vq  -  minimum  velocity  required  to  entrain  air  (1.1  m/sec). 
Elsawy  and  McKeogh  (1977)  found  for  circular  jets: 


1.2 


^3.7 


Fr, 


0.5 TTf*  ^ 


d®  FrJ 


1.66 


(D8) 


where 

V,  -  volume  of  entrained  air 
Frj  -  Froude  number  of  the  jet 

Fr^  -  Froude  number  of  the  jet  as  it  exits  the  nozzle 
d  -  diameter  of  the  jet 

4.  Several  slight  modifications  were  made  to  the  Elsway  equations  to 
make  them  applicable  to  weirs.  Combining  Equations  D7  and  D8,  an  equation  for 
the  deficit  ratio  r  can  be  written  as: 


r  =  exp 


(D9) 


where  k^a  is  the  liquid  film  coefficient.  From  this  Thene  developed  a  rela¬ 
tion  for  the  deficit  ratio  of: 


In  r  =0.156  Frf  ®®  ^ 


(DIO) 


which  can  be  transformed  into  an  equation  for  transfer  efficiency  at  20  “C  of: 


=■20 


1  -  exp  -  p . 156  Frt 


,2.69 


t^ 


(Dll) 


D9 


where  Fr^  is  the  Froude  number  of  the  flow  at  Impact  into  the  tailwater. 
Figures  D8-D11  describe  the  performance  of  the  equation  of  Thene  (1988). 

Nakasone  Equations 

5.  Nakasone  (1987)  attempted  to  combine  the  variables  of  fall  height, 
discharge,  and  tailwater  depth  to  characterize  the  transfer  efficiency  at 
weirs.  By  quantifying  the  effects  of  each  of  these  parameters  on  aeration 
efficiency,  Nakasone  identified  zones  where  reaeration  could  be  predicted,  and 
derived  the  set  of  four  equations. 


a  .  For  (D  +  l-5Hg)  <  1.2  m  and  q  235m®/h/m 
In  r2o  =  0.0785  (D  +  1 . 


(D12) 


where 

D  -  drop  height,  m 

Ho  -  critical  water  depth  on  the  weir,  m 
q  -  discharge  per  unit  width  of  weir,  m®/h/m 
r2o  “  deficit  ratio  at  20  ®C 
H  -  tailwater  depth,  m 
or 


E20  =  1  -  exp  -[0.0785  (D  ♦  1 . 

b.  For  (D  +  1.5Ho)  >  1.2  m  and  q  <  235m® /h/m 

In  r2o  =  0.0861  (D  +  1 . 5Ho)°  (D14) 

or 

E20  =  1  -  exp -[0. 0861  (D  ♦  1.5Ho)°  ®-^®q°  *^®H0  ®i°] 
or 


DIO 


Predicted  E 


Figure  D8.  Thene  equation  for  ogee  spillways. 
Predicted  values  of  gas  transfer  versus  measured  values 


Measured  E 

Figure  D9.  Thene  equation  for  gated  conduits. 
Predicted  values  of  gas  transfer  versus  measured  values 


Dll 


Figure  DIO,  Thene  equation  for  weirs. 
Predicted  values  of  gas  transfer  versus 


measured  values 


Measured  E 

Figure  Dll.  Thene  equation  for  gated  sills. 
Predicted  values  of  gas  transfer  versus 
measured  values 


D12 


c.  For  (D  -t- l.SHg)  S  1.2m  and  q  >  235m’/h/n 
In  r^o  -  5.39  (D  ♦ 

Ejo  -  1  -  exp -[5.39  (D  ♦ 

d.  For  (D  +  I.5H5)  >  1.2  m  and  q  >  235m®/h/m 

In  rjo  -  5.92  (D  ♦  i.5He)0  ««q-o.363H0.3io 
or 

E20  -  1  -  exp -[5. 92  (D  ♦  1.5Hc)'’  "^“q‘“  ®®®H°  ®^°] 


(D16) 

(D17) 


(D18) 


(D19) 


6.  The  results  of  these  equations  are  described  in  Figures  D12-D14. 

Holler  Equation 

7.  Holler  (1970)  used  a  form  of  Equation  7,  main  text,  and  proposed 
that  the  area  of  the  air/water  interface  is  a  function  of  momentum  change 
occurring  in  a  hydraulic  Jump  or  in  a  jet  impinging  on  a  water  surface.  Using 
a  series  expansion  for  the  exponential  term  in  Equation  6  and  reducing  the 
variables,  Holler  approximated  the  deficit  ratio  by 

r  »  1  ♦  ^(Av)2 

Av  is  then  related  to  the  head  by  energy  equilibrium,  and  for  weirs  and  over- 
falls  Holler  presented  the  equation  in  the  following  form: 


r  -  ^  +  1 


(D21) 


D13 


Predicted  E 


Measured  E 

Figure  D12,  Nakasone  equation  for  ogee  spillways. 
Predicted  values  of  gas  transfer  versus  measured  values 


Measured  E 

Figure  D13.  Nakasone  equation  for  weirs. 
Predicted  values  of  gas  transfer  versus  measured  values 
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Figure  D14.  Nakasone  equation  for  gated  sills. 

Predicted  values  of  gas  transfer  versus 
measured  values 

where  /3  was  equal  to  0.065  ft"^.  This  can  be  transformed  into  an  equation 
for  transfer  efficiency  at  20  *C: 


0.21325H 
(0.21325H  +  l) 


(D22) 


where  the  head  is  described  in  metres.  The  results  of  the  equation  of  Holler 
for  the  field  data  are  presented  in  Figures  D15-D18. 


8.  Tslvoglou  and  Wallace  (1972)  proposed  a  simple  mathematical  model 
that  predicts  gas  transfer  as  a  function  of  total  energy  dissipation.  For 
streams  the  model  takes  the  following  general  form: 

K,  -c^  (D23) 
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Figure  D15.  Holler  equation  for  ogee  spillways. 
Predicted  values  of  gas  transfer  versus  measured  values 


Measured  E 

Figure  D16.  Holler  equation  for  gated  conduits. 
Predicted  values  of  gas  transfer  versus  measured  values 
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Figure  D17.  Holler  equation  for  weirs. 
Predicted  values  of  gas  transfer  versus 
measured  values 


Measured  E 

Figure  D18.  Holler  equation  for  gated  sills 
Predicted  values  of  gas  transfer  versus 
measured  values 


where 


K2  -  reaeratlon  coefficient 
c  -  escape  coefficient  per  foot  of  energy  loss 
AE  •  energy  difference  between  upstream  and  downstream  points 
tf  -  time  of  flow  between  the  two  points. 

Combining  this  equation  with  the  first-order  reaeratlon  equation  of  Streeter 
and  Phelps  (1925) ,  Tslvoglou  and  Neal  (1976)  presented  the  equation  In  the 
following  form; 

E20  -  1  -  exp  (-^h)  (D24) 

where  ^  equals  0.054  ft~^  and  h  Is  the  head  In  feet,  or: 


E20  =  1  -  exp  {-0.1772h) 


(D25) 


where  h  Is  described  In  metres.  The  results  of  the  Tslvoglou  and  Wallace 
equation  are  presented  In  Figures  D19-D22. 

Foree  Equation 

9.  The  equation  of  Foree  (1976)  Is  Intended  to  characterize  the  aera¬ 
tion  that  occurs  over  small  falls  and  dams  found  In  creeks.  Foree  took 
Tslvoglou  and  Neal's  relationship  and  changed  the  coefficient  p  to  fit  his 
data  from  hydraulic  structures.  The  efficiency  at  20  ®C  can  be  written  as: 

E20  -  1  -  [exp  (-0. 5249h)f  ^^^6) 

with  head  measured  in  metres.  The  applicability  of  the  Foree  equation  to  low- 
head  structures  on  rivers  is  described  by  Figures  D23-D26. 

Preul  and  Holler  Equation 

10.  The  reaeratlon  prediction  equation  developed  by  Preul  and  Holler 
(1969)  was  Intended  to  be  a  generic  gas  transfer  equation  for  low  run  of  the 
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Figure  D19.  Tslvoglou  and  Wallace  equation  for  ogee  spillways 
Predicted  values  of  gas  transfer  versus  measured  values 
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Figure  D20.  Tslvoglou  and  Wallace  equation  for  gated  conduits 
Predicted  values  of  gas  transfer  versus  measured  values 
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Figure  D21.  Tsivoglou  and  Wallace  equation  for 
weirs .  Predicted  values  of  gas  transfer  versus 
measured  values 
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Figure  D22.  Tsivoglou  and  Wallace  equation  for 
gated  sills.  Predicted  values  of  gas  transfer 
versus  measured  values 
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Figure  D23.  Foree  equation  for  ogee  spillways. 
Predicted  values  of  gas  transfer  versus  measured  values 
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Figure  D24.  Foree  equation  for  gated  conduits. 
Predicted  values  of  gas  transfer  versus  measured  values 
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Figure  D25.  Foree  equation  for  weirs. 
Predicted  values  of  gas  transfer 
versus  measured  values 
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Figure  D26.  Foree  equation  for  gated  sills. 
Predicted  values  of  gas  transfer  versus 
measured  values 


D22 


river  dams.  Among  the  many  parameters  that  they  identified,  Preul  and  Holler 
narrowed  the  scope  of  Interest  to  Just  five  terms,  which  in  turn  could  be 
related  to  Just  one  term.  By  using  the  Froude  number  of  the  flow  before  the 
hydraulic  Jump,  a  relation  for  the  deficit  ratio  is  presented  as: 

rjo  -  1  +  666  ”  <>^27) 

where  is  the  Froude  number  at  Impact.  This  is  rewritten  in  terms  of 

efficiency  at  20  "C  as: 


-20 


jl  ♦  666 


(D28) 


Figures  D27-D30  present  the  comparison  of  the  predicted  values  to  actual 
measured  values  of  transfer  efficiency. 


Rindels  and  Gulliver  Equation 


11.  Rindels  and  Gulliver  (1991)  presented  an  equation  of  the  form: 


where  the  values  are  all  in  non-SI  units,  or: 


E20 


1  -  exp 


f  -0.2625h 

[1+  0.2153q 


(D29) 


(D30) 


where  Zp  is  equal  to  the  tailwater  depth,  h  is  the  head  across  the 
structure,  and  q  is  the  discharge  per  unit  width  (all  values  in  SI  units). 
The  intention  was  to  further  refine  the  Tsivoglou  and  Wallace  (1972)  equation 
by  incorporating  the  effects  of  tailwater  pool  depth  on  aeration  with  the 
effects  of  spillway  height.  The  authors  reasoned  that  transfer  efficiency 
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Figure  D27.  Preul  and  Holler  equation  for  ogee  spillways 
Predicted  values  of  gas  transfer  versus  measured  values 
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Figure  D28.  Preul  and  Holler  eqxiation  for  gated  conduits 
Predicted  values  of  gas  transfer  versus  measured  values 
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Figure  D29 .  Preul  and  Holler  equation  for  v’airs. 
Predicted  values  of  gas  transfer  versus  measured 

values 


Measured  E 

Figure  D30.  Preul  and  Holler  equation  for  gated 
sills.  Predicted  values  of  gas  transfer  versus 
measured  values 


would  be  proportional  to  the  depth  of  tailwater  and  height  of  spillway  and 
would  be  Indirectly  proportional  to  the  unit  discharge. 

12.  The  results  of  the  equation  of  Rlndels  and  Gulliver  are  presented 
In  Figures  D31-D33. 


13.  Wilhelms  examined  the  hydraulic  conditions  at  gated  sill  structures 
and,  using  an  exponential  form  of  equation,  empirically  related  the  oxygen 
transfer  to  the  discharge,  head  loss  across  the  structure,  and  submergence  of 
the  gate  lip.  The  exponent  In  the  gas  transfer  equation  was  directly  related 
to  the  discharge  and  head  loss  and  Inversely  related  to  submergence: 

E  20  -  1  -  exp  I -0.000797  ^  -  0.188  J  ^^^1) 

where 

h  -  head  across  the  structure,  ft 
q  -  unit  discharge,  ftVsec 
s  -  submergence  of  gate  lip,  ft 

Caution  was  advised  In  using  the  equation  because  of  prediction  errors  for 
very  small  submergences  or  for  large  discharges  that  change  the  hydraulic 
action  In  the  structure  stilling  basin.  Figures  D34-D36  show  the  results  for 
the  Wilhelms  equation. 

Wilhelms  and  Smith  Adjustment  to  Tslvoglou  Equation 

14.  Wilhelms  and  Smith  (1981)  modified  the  coefficient  in  the  equation 
of  Tslvoglou  and  Wallace  (1972)  for  a  data  set  that  consists  mostly  of  gated 
conduits.  For  their  data  set  the  following  relation  best  predicted  the 
transfer  efficiency: 

E  =  1  -  exp  (-0.  1476h) 


where  h  Is  the  head  across  the  structure.  Figures  D37-D40  present  the 
results  of  the  Wilhelms  and  Smith  equation  for  several  types  of  structures. 
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Figure  D31.  Rlndels  and  Gulliver  equation  for  ogee  spillways. 
Predicted  values  of  gas  transfer  versus  measured  values 
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Figure  D32.  Rindels  and  Gulliver  equation  for 
weirs.  Predicted  values  of  gas  transfer  versus 
measured  values 
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Figure  D33.  Rindels  and  Gulliver  equation  for 
gated  sills.  Predicted  values  of  gas  transfer 
versus  measured  values 
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Figure  D34.  Wilhelms  equation  for  ogee  spillways. 
Predicted  values  of  gas  transfer  versus  measured  values 
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Figure  D35.  Wilhelms  equation  for  weirs. 
Predicted  values  of  gas  transfer  versus 
measured  values 
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Figure  D36.  Wilhelms  equation  for  gated  sills. 
Predicted  values  of  gas  transfer  versus  measured 

values 
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Figure  D37.  Wilhelms  and  Smith  equation  for  ogee  spillways. 
Predicted  values  of  gas  transfer  versus  measured  values 


Measured  E 

Figure  038.  Wilhelms  and  Smith  equation  for  gated  conduits. 
Predicted  values  of  gas  transfer  versus  measured  values 
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Figure  D39.  Vilhelms  and  Smith  equation  for 
weirs.  Predicted  values  of  gas  transfer 
versus  measured  values 


Measured  E 

Figure  040.  Vilhelms  and  Smith  equation  for 
gated  sills.  Predicted  values  of  gas  trans¬ 
fer  versus  measured  values 
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APPENDIX  E:  NOTATION 


El 


a  Specific  surface  area,  A/V 

A  Surface  area  associated  with  the  volume  V  ,  over  which  transfer 
occurs 

Bg  Bias  uncertainty  in  the  measurement  of  and  Cf  ;  bias 

uncertainty  In  DO  concentration  due  to  the  tltrant 

Be  Bias  uncertainty  In  C, 

Bg  Bias  uncertainty  in  the  transfer  efficiency 
C  Ambient  concentration  in  water 

Ccj^  Concentration  of  chloride  ions,  g/i 

Cj  Final  oxygen  concentration 
Ci  Initial  oxygen  concentration 

C,  Saturation  concentration 

Cji  saturation  concentration  of  gas  i  ,  g/m^ 

dC  Change  in  concentration  across  the  interface  thickness 
dC/dt  Rate  of  change  of  concer  ^ration 

dD/dt  Rate  of  change  of  the  saturation  deficit 

ds  Interface  thickness 

D  Saturation  deficit 

Df,Di  Final  and  initial  oxygen  deficits,  respectively 

Dj,  Molecular  diffusion  coefficient 

E  Transfer  efficiency 

E  Best  estimate  of  the  true  value  of  transfer  efficiency  that  is 
predicted  by  the  equation 

E„,Ep  Measured  and  predicted  oxygen  transfer  efficiency,  respectively 

E20.  Ej  Gas  transfer  efficiencies  at  20  °C  and  at  T  °C 

Fj  Froude  number  of  the  jet 

g  Acceleration  due  to  gravity,  m/sec^ 

h  Head  loss,  m 

H  Tailwater  depth,  m 

Hg  Critical  water  depth  on  the  weir,  m 

Henry's  law  constant  for  gas  i  ,  g/m®  -  atm 
J  Mass  flux  per  unit  area  across  an  interface  or  a  plane 
kj.  Liquid  film  coefficient  for  oxygen  transfer 
K2  Reaeration  or  oxygen  transfer  coefficient 

Elevation  above  mean  sea  level,  m 


M 


Nf 

P 

^atm 

Pi 

q 

Q 

r 

R 

Rh 

s 

t 

t* 


T 

U 

Ue 

V 


Wc^.Wc^ 

We 

1/ 


b,p,t,v 


d 

D 

E 

E^ 

fj 

Fr^ 


Nu’.aber  of  upstream  and  downstream  measurements,  respectively 
Froude  number  at  the  point  of  impact 
Probability  of  the  confidence  interval  on  efficiency 
Reported  pressure,  mm  of  mercury 

Partial  pressure  of  gas  1  in  the  atmosphere,  atm 

Specific  discharge,  volume/time/length;  discharge  over  the 
structure 

Total  discharge  through  a  structure 
Deficit  ratio 
Reynolds  number 

Hydraulic  radius  of  the  flow  in  the  conduit 
Gate  submergence ,  m 

Elapsed  time  from  initial  to  final  deficits 

Student's  t-value  corresponding  to  nc  -  1  or  nc  -  1  degrees  of 
freedom  and  a  95  percent  confidence  interval 

Temperature 

Uncertainty 

Total  uncertainty  in  the  transfer  efficiency 
Volume  of  the  water  body  over  which  A  is  measured 
Precision  uncertainties  in  and  ,  respectively 
Precision  uncertainty  in  the  transfer  efficiency 
Kinematic  viscosity 

Standard  deviation  of  the  upstream  and  downstream  measurements, 
respectively 


Appendix  D 

Width,  perimeter,  thickness,  and  velocity  of  the  jet, 
respectively 

Escape  coefficient  per  foot  of  energy  loss 
Diameter  of  the  jet 
Drop  height,  m 
Efficiency 

Gas  transfer  efficiency  at  a  given  temperature  T 
Indexing  coefficient 

Froude  number  of  the  flow  at  Impact  into  the  tailwater 
Froude  number  of  the  jet 


EA 


Fr^  Froude  number  of  the  jet  as  it  exits  the  nozzle 
g  Acceleration  due  to  gravity 
h  Head  across  the  structure;  fall  hel^t 
H  Depth  of  the  tallwater  pool 
Hq  Critical  water  depth  on  the  weir 

K^a  Liquid  film  coefficient 

K2  Reaeratlon  coefficient 
Nf  Froude  number  at  Impact 

q  Discharge  over  the  structure;  discharge  per  unit  width  of  weir, 
m^/h/m 

Q,  Air  discharge 

Q„  Water  discharge 

r  Deficit  ratio 

ri5  Deficit  ratio  at  15  ®C 

r2o  Deficit  ratio  at  20  *C 

R  Reynolds  number  of  the  jet 
s  Submergence  of  gate  lip,  ft 
tj  Time  of  flow  between  the  two  points 
Vg  Minimum  velocity  required  to  entrain  air  (1.1  m/sec) 

Vg  Volume  of  entrained  air 

Zp  Tallwater  depth 

AE  Energy  difference  between  upstream  and  downstream  points 
1/  Kinematic  viscosity 
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